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Abstract 

The polymer systems are discussed in the framework of the Landau-Ginzburg 



O ■ 

model. The model is derived from the mesoscopic Edwards hamiltonian via the 



conditional partition function. We discuss flexible, semiflexible and rigid polymers. 
The following systems are studied: polymer blends, flexible diblock and multi-block 
copolymer melts, random copolymer melts, ring polymers, rigid-flexible diblock 
copolymer melts, mixtures of copolymers and homopolymers and mixtures of liquid 
crystalline polymers. Three methods are used to study the systems: mean-field 
model, self consistent one-loop approximation and self consistent field theory. The 
following problems are studied and discussed: the phase diagrams, scattering inten- 
sities and correlation functions, single chain statistics and behavior of single chain 
close to critical points, fluctuations induced shift of phase boundaries. In particular 
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we shall discuss shrinking of the polymer chains close to the critical point in poly- 
mer blends, size of the Ginzburg region in polymer blends and shift of the critical 
temperature. In the rigid flexible diblock copolymers we shall discuss the density 
nematic order parameter correlation function. The correlation functions in this sys- 
tem are found to oscillate with the characteristic period equal to the length of the 
rigid part of the diblock copolymer. The density and nematic order parameter mea- 
sured along the given direction are anticorrelated. In the flexible diblock copolymer 
system we shall discuss various phases including the double diamond and gyroid 
structures. The single chain statistics in the disordered phase of flexible diblock 
copolymer system is shown to deviate from the Gaussian statistics due to fluctua- 
tions. In the one loop approximation one shows that the diblock copolymer chain 
is stretched in the point where two incompatible blocks meet but also that each 
block shrinks close to the microphase separation transition. The stretching out- 
weights shrinking and the net result is the increase of the radius of gyration about 
the Gaussian value. Certain properties of the homopolymer, copolymer system are 
discussed. The diblock copolymers solubilizes two incompatible homopolymers by 
forming a monolayer interfaces between them. The interface has the positive sad- 
dle splay modulus which means that the interfaces in the disordered phase should 
be characterized by the negative Gaussian curvature. We also show that in such 
a mixture the Lifshitz tricritical point is encountered. The properties of this un- 
usual point are presented. The Lifshitz, equimaxima and disorder line are shown 
to provide a useful tool for studying local ordering in polymer mixtures. In the 
liquid crystalline mixtures the isotropic nematic phase transition is discussed. We 
concentrate on static, equilibrium properties of the polymer systems. 
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I. Introduction 

The introduction of the order parameter in physics was a milestone in the 
description of phase transitions. The same Landau-Ginzburg approach has been 
applied to such diverse phenomena as transitions in liquid crystals and superconduc- 
tivity, liquid-solid or solid-solid phase transitions and transitions in the superfluid 
helium, ferromagnetic paramagnetic phase transition and phase transitions in poly- 
mer systems. The purpose of this review is the application of the Landau-Ginzburg 
model to polymer systems. Despite the obvious limitations of this approach (valid 
for weakly first order and continuous phase transitions) it provides plenty of infor- 
mation about the polymer systems such as: polymer blends, diblock copolymers 
and liquid crystalline polymers. In polymer systems the parameters characterizing 
the Landau-Ginzburg free energy can be directly related to the mesoscopic parame- 
ters of the polymer systems. The mesoscopic parameters characterizing for example 
the interactions are, on the other hand, related to the microscopic interaction pa- 
rameters via the integrals of the direct corelation functions. Therefore the results 
obtained in the Landau-Ginzburg model shed light on the relation between the 
mesoscopic and/or microscopic parameters of the system and its behaviour. 

The simplest polymer mixture is the A homopolymer, B homopolymer blend. 
The very fruitful approach to the study of the system near its consolute (critical) 
point is based on the mean field Flory-Huggins free energy 1 ) and the de Gennes 
random phase approximation (RPA) 2 ) for the scattering intensity. Both can be 
stated in terms of the Landau-Ginzburg model as we show in Section III. RPA and 
Flory-Huggins approach are two compatible mean-field models. They are named 
mean-field since in both models only the most probable configurations of the poly- 
mer chain are taken into account and, more importantly, the critical long-wavelength 
fluctuations are neglected. It was de Gennes who pointed out 3 - 1 that the mean-field 
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theory is rather good for high molecular mass polymer mixtures, in contrast to the 
low molecular mass mixtures, for which the mean-field theory breaks down close 
to the critical point. In fact, in the limit of Na,Nb — ► oo, the mean-field theory 
is exact 3-5 ). The mean-field approach is quantitatively correct if the correlation 
length £ ~ (T/T c — l)" 1 is comparable or smaller than the typical length scale in 
the system i.e. the radius of gyration ~ \fN . Comparing two length scales, for 
large N, we find that the mean-field theory breaks down very close to the critical 
point, i.e. for 

|T - T c \/T c ~ 1/N. (1.1) 

This is the Ginzburg criterion. The region around the critical point where the mean 
field theory breaks down is called the Ginzburg region. The shift of the critical 
temperature, T c , with respect to T™f (calculated in the Flory Huggins model) is 

\T C - T™ f \/T™ f ~ 1/y/N (1.2). 

The shift, induced by fluctuations, is obtained in the one-loop approximation 6 ' 7 ''. 
It is larger than the Ginzburg region itself. The derivation of Eq(1.2) is based on 
the study of the scattering intensity, I(q), where q is the scattering wavevector. 
For T -> T c and q -> 0, we expect I _1 (0) ~ (T/T c - 1)-t, where 7 = 1 in the 
mean- field theory 6 ' 7 ). The breakdown of this relation is usually used in experiment 
to estimate the Ginzburg region 7 ' 8- 13 ) and therefore, we also base the calculations 
of the Ginzburg region on the scattering intensity 14 ' 15 ). Sufficiently close to T c , 
fluctuations become important and 7 changes from 1 to its universal value of 1.26. 
This change can be extracted if one plots log 1(0) versus log \T/T C — 1|. 

Different polymers can be combined into single material in many ways and 
polymer blends find widely spread application. Many of the potential application 
depends on the morphology of the system which in turn depends on the single 
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chain behaviour. Here we shall discuss the behavior of the single chain in the 
polymer blend 16-19 ) as well as in the diblock copolymer system 20 ' 21 ). The chains 
shrink considerably at the approach to the critical point in homopolymer blends. 
Intuitively we can expect that when the correlation length becomes larger than the 
radius of gyration the chains demix on this scale 18 ) . On the other hand the chains 
stretch in diblock copolymer system 20 ' 21 ). 

An A-B diblock copolymer is a polymer consisting of a sequence of A-type 
monomers chemically joined to a sequence of B-type monomers. Even a small 
amount of incompatibility (difference in interactions) between monomers A and 
monomers B can induce phase transitions in a mixture of homopolymers. However 
A-homopolymer and B-homopolymer are chemically joined in a diblock, so that 
a system of diblocks cannot undergo a macroscopic phase separation. Instead a 
number of order- disorder phase transitions take place in the system between the 
isotropic phase and spatially ordered phases in which A-rich and B-rich domains, 
of the size of a diblock copolymer, are periodically arranged in lamellar, hexago- 
nal, bcc, lamellar-catenoid, double diamond (Fig.l) or gyroid structures (Fig. 2). A 
Landau- Ginzburg model of the phase transitions in diblock copolymer system was 
first formulated by Leibler 22 ) and later refined by Fredrickson et al 20 ' 23 ). Matsen and 
Schick 24 ' 25 ) applied the self consistent field theory 26 ' 27 ) to the phase transitions in 
this system in the weak segregation limit 28 ) , where the domain boundaries between 
the A rich and B rich domains are wide. The opposite case is the strong segregation 
regime with very sharp and narrow domain boundaries; here we consider only the 
former case. The theoretical studies of the system were greatly stimulated by exper- 
iments. The cylindrical and spherical and lamellar phases have been known for long 
time 29 ), but only recently the novel bicontinuous structures have been discovered 
in the diblock copolymer systems 30-32 ) (Fig.l). In the polystyrene-polyisoprene 
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system (Fig. 3) the bicontinuous double diamond structure consists of two channels 
each of diamond symmetry separated by the surface (Fig.l), which as has been 
argued by Thomas et al 32 ) is of constant mean curvature 33 ) . 

The richness of the phase diagram and behavior of the polymer system is en- 
hanced if we mix the diblock copolymers and homopolymers 34-37 ). Ternary mix- 
tures of A,B homopolymers and A-B diblock copolymer 34-36 ) are similar in some 
respects to mixtures of oil, water and amphiphile 38 ) . The copolymer acts as a sur- 
factant for the A,B homopolymer blend, decreasing the surface tension between 
A-rich and B-rich phases as it accumulates on the interface between these two 
phases and reduces the number of energetically unfavorable contacts between the 
A-homopolymers and B-homopolymers. In spite of this similarity, the comopolymer 
is not generally as efficient a solubilizer of homopolymers as a good amphiphile is 
of oil and water 35 ). The analysis of the disordered system is based on the study 
of disorder 39 ), Lifshitz and equimaxima lines. The disorder line of the system is 
the locus of points at which all correlation functions in general no longer decay 
monotonically, but contain an exponentially damped oscillatory component. This 
oscillation reflects the tendency of the copolymer to order the A (B) monomers of 
the system. The Lifshitz line is the locus of points at which the peak in the structure 
function just begins to move off of zero wavevector. It therefore indicates the point 
at which oscillatory components begin to dominate the particular correlation func- 
tion, in contrast to the disorder line which indicates the point at which oscillatory 
components appear in all correlation functions. The Lifshitz line of the structure 
function of all A monomers is quite close to the disorder line, indicating that once 
the tendency to order the system appears, it quickly dominates the response of this 
function. In contrast, the Lifshitz line associated with the structure function of 
those monomers located only on the homopolymers is, in a large part of the phase 
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diagram, far from the disorder line. These results taken together indicates that the 
copolymer is rather ineffective in organizing the homopolymers, and much more ef- 
ficient in organizing itself, an organization expressed most strongly in the lyotropic 
phases which easily exist in the complete absence of homopolymer. Beyond the 
disorder line the interface between the A-rich homopolymer B-rich homopolymer 
phase is not wetted by the disordered phase 36 ) so it consists of of a thin copolymer 
layer. In the disordered phase we expect that A-rich and B-rich domains are also 
separated by the thin copolymer layer. Assuming that the structure of the interface 
is the same as the internal interfaces in the disordered phase Matsen and Schick have 
calculated 36 ) the profiles and from them obtained the saddle-splay modulus for the 
internal intefaces. The modulus is positive, thus internal interfaces in the disordered 
phase of A,B homopolymers and AB diblock copolymers are characterized by the 
negative Gaussian curvature. 

Random copolymers 40 -* are produced by polimerization reaction in the medium 
that contains two (or more) distinct monomer species. The chemical composition 
of the copolymers is dictated by a set of bimolecular reaction rates. The system 
is characterized by quenched disorder 41-43 ). The chain architecture is determined 
statistically by the conditional probability that A monomer follows B monomer. 
Under steady state conditions of the reaction this probability is independent of 
the location of the monomer in the chain. The structures formed in the random 
copolymer system depend on the structure of the chain and thus on the polimer- 
ization process, thus the phase diagram can be a useful guide for the design of new 
polymeric materials. 

In the multiblock copolymer system 44 ) new ingredients enter into physics in 
comparison to diblock copolymers. In the latter the junction of the blocks of a 
given copolymer is constrained to one A/B block interface; in the former bridging 
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and looping is possible. In the loop configuration of the AB multi-block copolymer 
system both junctions can reside on the same interface (Fig. 4). In the bridge con- 
figuration the middle B block spans over one B-rich lamella and joins two A-rich 
lamellas. In many material properties the bridge configurations are important 45 ) . 
The bridging fraction thus provides important information for the design of new 
materials. 

Ring homopolymers or ring diblock copolymers provide yet another example of 
the influence of chain architecture on the structure and phase diagram of polymer 
system. In the AB ring polymer the A block is connected with the B block at two 
junction points, therefore the chain has no free ends. The topological constraint 
introduced by this chain architecture, prevents the different chains from entangled- 
ment or self-knotting (Fig. 5). Thus the statistics of the rings is not Gaussian and 
the size of a ring is intermediate between the Gaussian and collapsed size. For the 
ring homopolymer in the melt the scaling arguments 46 "-* show that the radius of 
gyration scales with the polimerization index, AT, as N u , with v = 0.4. The appli- 
cation of RPA to the Edwards hamiltonian with the specific constraint of no knots 
gives the more accurate value 46b ) of v = 0.445. For the ring diblock copolymer the 
computer simulations 47 ) give v = 0.45. Please note that at high temperature the 
scaling in the melt of ring polymers should be the same as in the ring diblock copoly- 
mer system. The random phase approximation, with assumed Gaussian statistics, 
shows 48 ' 49 ) that the transition to the ordered lamellar phase occurs at the temper- 
ature T ring = 1.78Tu near (for symmetric rings), where Ti inear is the temperature 
for the transition to the lamellar phase in the linear diblock copolymer system. 
This result compares very well with simulations 47 ), despite the fact that the ring 
statistics is certainly not Gaussian. When we decrease the temperature the rings 
strongly stretch along the axis connecting the centers of mass of two blocks. 
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So far we have discussed the flexible-flexible diblock copolymer systems, where 
the diblock copolymers are modelled as Gaussian chains or, almost equivalently, 
chains with freely rotating bonds of fixed length. One can investigate the changes 
in the structure and stability of the disordered (isotropic) phase of the diblock 
copolymer system with the stiffness of a copolymer molecule. Microphase separated 
block copolymers with blocks of different stiffness are technologically important as 
composites. The flexible part of the polymer provides resistence to fracture while 
the rigid one resistence to tensile stresses and thermal stability. The simplest system 
of rigid-flexible diblock copolymers is discussed here 50 -*. In a molecule one part of 
the copolymer is a chain with freely rotating bonds while the other is the rigid 
rod. In such a system, the density of monomers alone is not enough to specify the 
structure of the rod-rich domains nor the stability limits. In order to have a complete 
description of the system we must include a nematic order parameter tensor 51 ). In 
the high temperature isotropic phase, the various correlation functions between 
density and nematic order operators provide information about the structure. The 
correlation functions, in the Gaussian approximation, are found to oscillate with 
a characteristic period equal to the length of the rigid part. The lamellar domain 
boundaries of systems with longer rods are sharper than those with shorter ones, 
and the linear density measured along any direction is anticorrelated with the order 
parameter measuring nematic order in that direction. This latter effect is due to 
the isotropy of the system, and is relatively independent of the Maier-Saupe and 
Flory-Huggins parameters. This model is a special example of the more general 
model 52 ) shown in Fig. 6. Here the A block consists of segments of small size 6^ and 
B block consisting of segments of much larger size bs- 

Liquid crystal polymers and their mixtures are studied both because of their 
practical utility and their intrinsic interest 53 ' 54 -*. In the case of mixtures 55 ' 56 ), one 
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wants to know how the location of various phases, isotropic and nematic, and their 
transitions depend on the properties of the two components, their rigidities, poly- 
merization indices, interactions etc. An approach to this problem requires a model 
for the liquid crystal polymers. Recently, Maissa and Sixou 57 ) addressed these is- 
sues on by modelling the polymers as elastic lines with an energy of bending 58 ). 
Their calculation can be applied to many different mixtures, and they obtained 
various phase diagrams numerically for several values of the interaction constants. 
It is difficult, however, to extract general relationships between properties of the 
phase diagrams and the polymer parameters from their formalism. Such depen- 
dences can be made manifest if the Landau de Gennes expansion of the free energy 
is derived from the polymer model because the coefficients in this expansion depend 
upon the polymer characteristics of interest. Fredrickson and Leibler 59 ), inspired 
by experiments of Moore and Stupp 56 ) , carried out this program in a model which 
also treated both polymers as elastic lines, ones characterized by different bending 
constants. Among other results, they found that the expansion was only applicable 
to this model if the two bending rigidities were large and nearly equal. On the other 
extreme there is a model of rigid flexible mixture where the differences in stiffness 
are large 60 -*. 

There are many phenomena which will not be discussed in this article. The 
mixtures are assumed as incompressible and compressiblity effects are ignored 61 ). 
The crossed-linked blends 62 ) and charged polymer systems 63 ) , although interesting, 
will not be discussed here. We have concentrated in this paper on the equilibrium 
properties of various polymer systems; consequently the dynamics of blends 64-66 ) 
and diblock copolymer system 67 ' 68 ) is not included. 

The paper is organized as follows: In section II the mesoscopic model of poly- 
mer molecules is presented. Here we discuss the bond configuration for different 
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systems, the various interaction hamiltonians steming from the Edwards model, 
density operators and order parameters and finally the conditional partition func- 
tion. In section III we present the detailed study of the polymer blend system, 
including the discussion of the critical temperature in the Gaussian and one-loop 
approximation, the single chain statistics and the problem of the upper wavevector 
cutoff in polymer blends and in general in polymer systems. In section IV we discuss 
the diblock copolymer system. Here we discuss the behaviour of tagged chain in the 
melt of copolymers, phase diagram etc. The system is studied in the random phase 
approximation, one-loop approximation and self consistent field theory. In section 

V we study the random copolymer and multiblock copolymer systems. In section 

VI we present the results for the mixtures of homopolymers and copolymers. The 
correlations in the rigid-flexible diblock copolymer system are studied in section VII 
in the random phase approximation. The mixtures of rigid and flexible polymers 
are discussed in section VIII. 

II. The mesoscopic model of a polymer system 
The mesoscopic model of polymers is specified by the Hamiltonian describing 
the interaction between monomers, and the density distribution of monomers within 
a polymer molecule. Furthermore one has to specify the density operators which 
are needed to specify the structure of the system. The average of this operators is 
related to the order parameters. Next the conditional partition function has to be 
constructed. This is the partition function for the system subject to the constraint 
of the the fixed configuration of the order parameter field. The logarithm of this 
partition function expanded in the order parameter gives the Landau- Ginzburg free 
energy. The calculations of correlations are also based on this partition function. 
Below we specify all the elements of the model. 
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A. Bond configuration 

The information that N monomers are connected and form a chain is specified 
by the distribution function W[r]. To describe the flexible polymer, we use the 
chain model, in which atoms are described as being joined by freely rotating bonds 
of fixed length /. The distribution for N atoms in such a chain is given by 69 ) 

N 



and is normalized as 



with 



DrW[r] = 1, (2.2) 



Dr=^dr dr 1 ---dr N . (2.3) 

Here Uj = r~ rj_i is a vector specifying the orientation of the monomer and is the 
location of the i point between the subsequent monomers. The same quantitative 
results at the level of the Landau- Ginzburg free energy are obtained if instead of 
the previous model we use the Gaussian model in the discrete version 

N 

21' 



WM-expj-A^lr.-r.^ j (2 . 4 ) 



or a continuous version of the model, 



N 



ds 



dr(s) 



ds 



2 N 



(2.5) 



In both cases one normalizes the distribution. 

The rigid polymers are described as rigid rods, or needles, in which all bonds 
have not only the same length, but also the same direction. In this case 

N 

47T/ 2 



^W = ^^ il n^- U ^)- (2-6) 

J=2 
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Finally, for a diblock copolymer in which a flexible chain consisting of N atoms 
has been joined with a rod containing M atoms, 

N o/i | _ n N+M 

^'1 = 11^^ II «(«,-«,-.) (") 

i=l j=iV+l 

For semiflexible polymers one can use the model shown on Fig. 6 which in the 
extreme limit reduces to the one discussed above. A different model of semiflexible 
polymers is the wormlike-chain model 58 ' 59 ). In the discrete version we have: 

'l N 



W[r] ~ exp l^-h u.+iu.J , (2.8) 

where k\ is the bending elastic constant, which is the energy penalty for the change 
of angle between the nearests monomers. The continuous version of this model 58 ) 
is defined as follows: 



W[r] ~ exp I ~2^ 1 I ^ s 



1, f N . du(s) 







2 X 



(2.9) 



ds 

subject to the constraint |u| 2 = 1. The next step in the definition of the model is 
the specification of the interactions between the monomers. 

B. Interactions between monomers 
The typical length scale for a flexible chain is given by the radius of gyration i.e. 
the size of the region occupied by the chain in the melt. In the simplest Gaussian 
approximation one finds that it is proportional to \/Nl and is much larger than 
the monomer size /. All the interesting phenomena take place at the length scale 
proportional to the radius of gyration. On the other hand the range of the potential 
is proportional to /, and therefore is not relevant to the phenomena discussed in 
the paper. Guided by this simple observation the following short range effective 
interaction potential has been proposed by Edwards 70 - 1 : 

v if( r i, r j) = Wa^fa - Tj), (2.10) 
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where 6 is the Dirac delta function and w a p is the effective interaction parameter 
for a and (3 type monomers (e.g. A, B monomers in the binary homopolymer 
mixture). The effective interaction parameter is given by the integral 70 ' 71 -* of the 
direct correlation function 72 ' 73 \ In the first approximation for the rigid molecules, 
such integral, does not depend on temperature and is equal to the excluded volume, 
the volume inaccesible to one molecule when the other is fixed in space. In the 
case of two spheres of radius D the excluded volume is the larger sphere of radius 
2D. Apart from the repulsive forces there is also an attractive potential e.g. van 
der Waals potential. The direct correlation function is, in the first approximation 
proportional to this potential. Thus w a p contains two contribution: one associated 
with the excluded volume and one with the attractive potential. Summarizing: the 
interactions given by Eq.(2.10) are related to the microscopic interaction potential 
via the direct correlation function. The delta function signifies the extreme short 
range character of the potential as measured in terms of the length given by the 
radius of gyration. 

The interaction potential for rigid, elongated molecules, depends not only on 
the positions but also on the orientations of molecules. In the case of two rods 51 ) 
with fixed orientations the excluded volume has roughly the shape of a rectangular 
box. If we expand the excluded volume in terms of the Legendre polynomials of the 
cosine of the relative angle of two rods we obtain the constant term, proportionl to 
w a p, plus the second term related to the second Legendre polynomial, P^. If higher 
order terms are neglected the anisotropic part of the potential is modelled as 50 ' 60 ); 

v$(r u m, vj, uj) = v aP 5{ ri - rj )P 2 (t^) (2.11) 

V \ u i\ \ u j \ J 

The total potential is given by the sum of Eq(2.10) and (2.11). Please note that 
in general the parameters w a p and v a p are not independent, since they follow from 
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the same direct correlation function. 

Finally for the dipolar monomers one can model the potential as: 

v?f (Ti, ju, Tj, uj) = v af3 S{ri -tj)! ) (2.12) 

V \ u t\ \ u j \ J 

The next step is the specification of the density operators needed for the defi- 
nition of the order parameters. 

C. Density operators 

For the sake of clarity we shall discuss this point using the example of the rigid- 
flexible diblock copolymer system 50 ^. The chain architecture is given by Eq(2.6). 
The system consists of n chains each consisting of Na A-type monomers joined to 
Nb, B-type monomers. The interactions between two monomers are described by 
the sum of Eq(2.10) and Eq(2.11). The density operators needed to specify the 
configuration of the system are as follows: 

1 n N A 

A n) ( r ) = -£X>(r-r7); (2-13) 

P0 7 =1 i=l 
n N A +N B 

4 n V) = f£ E 5 ( r - r D; ( 2 - 14 ) 

7=1 i=l ^ III / 

cv> = f ± T B *(«■ - '?) (f (u?> r u ;f" 7> ' - ¥) • < 2 - 16 ) 

P° 1 =i 1 =n a +i vz |u * ! z ' 

The first two operators, 0^ (r) and (j)^ (r) (superscript (n) stands for the number of 
chains in the summation) , represent the microscopic number fractions at point r of A 
and B monomers respectively; the next two tensor operators, ( r ) and Q^g (r) 
(cm, P = 1,2,3), represent the microscopic nematic tensor order parameters 51 -* at 
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point r for the A and B monomers respectively. These tensors are symmetric and of 
zero trace, thus each has five independent components. The total number of inde- 
pendent operators is, therefore, twelve. Here po is the average density of monomers 
in the system. 

The total interaction hamiltonian, H, is obtained by summing all the interac- 
tions between different monomers. Using the density operators it can be rewritten 
in the following simple form: 

H = k B Tp J dr(±«; AA (A n) ) 2 (r) + ^^ 

(2.17) 

where summation over repeated a and (3 indices is implied. 

There is no general recipe for the choice of the density operators. The hint as 
to the right choice is provided by the form of the interactions between monomers 
and the chain architecture. Here $y , QaP 7 = A, B form a minimal set of order 
parameter operators necessary to describe the system in the Landau approach. 

D. Conditional partition function 

The conditional partition function, Z[Pj\ is the partition function for the system 
subject to the constraint that the microscopic operators, P t be fixed at prescribed 
values, Pi(r) i.e. 

Z[Pi] = (l[s[Pi-Pi)y (2.18) 
where the average is calculated as follows: 

(••')= z~ 11 / Drl ■ • • W ^ eM~H/k B T) (2.19) 

° 7=1 

Here k B T is the Boltzman factor and Z is the canonical partition function. The 
Landau- Ginzburg (LG) free energy (in the mean- field approximation), 0[Pj], is 
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given by 

Q.[P i ] = -k B T\nZ[P i \. (2.20) 

The conditional partition function and the LG free enrgy are functionals of the order 
parameters P{. The partition function of the system is given by the summation of 
the conditional partional function over all possible configurations of the fileds Pf. 

Z = Il[ DP^}. (2.21) 

In the following sections we shall apply this method to polymer systems. 
III. Binary mixture of flexible homopolymers 

We consider a mixture of n A A-type polymers with N A monomers in each 
molecule, and n# B-type polymers with monomers in each molecule inside a 
volume V. The anisotropic interactions, v a p are neglected. For later convenience 
we introduce two fields, 

y A (r)=^ A \r)-0 (3.1) 

and 

vMr)=<^ B) (r)-(l-0). (3-2) 

These fields describe the excess of the fraction of A or B monomers at the point r 
over their average values in the system, = uaNa/^aNa + nsNs), and 1 — 0; 
they are zero in the homogeneous system. Finally, we assume that the system is 
incompressible, i.e. 

^>(r) + ^(r) = l, (3.3) 

so that 

tt A (r) = -*s(r) = tt(r). (3.4) 

For the Fourier transform of \P(q) = j <ir\I/(r)e irq , it implies *(q = 0) = 0. Since 
the system is incompressible, only one interaction parameter is needed, namely the 
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dimensionless Flory-Huggins parameter, x, waa + wbb — 2wab = — 2% < 0. In 
general x is positive so that it is energetically favorable for the system to separate 
into A-rich and B-rich phases. The first aim of the study is the location of the 
critical point for the mixture, below which the system separates into A-rich and 
B-rich phases. 

The field \1/ is for the incompressible blend the only order parameter and con- 
sequently both the partition function and the LG free energy are the functionals of 
Please note that \& vanishes in the homogeneous system, but is nonzero at the 
coexistance curve. We find, using the results of the cumulant expansion given in 
Appendix A and B, the following form of for \1/ independent on position: 

^KK^ivb))-*)* 2 

12 \ Nb (1-4,)3 N A <jr>) 

Linear terms, in have been omitted in Eq.(3.5) as they simply correspond to a 
shift, /lip, in the exchange chemical potential fx = (/ja/Na) — (hb/Nb), i-e. 

hi 6 ln(l-0) 1 1 , . . 

Here ha, \±b are the chemical potentials of A and B monomers. Please note that 
Eq(3.5) can be also obtained (as it should) directly from the Flory-Huggins free 
energy by expanding it in The critical point for the mixture corresponds to 
the vanishing of the first and second term in Eq(3.5), or in other words to the 
vanishing of the second and third derivative of the Flory-Huggins free energy with 
respect to the concentration. The location of the critical point in this mean-field 
approach (Appendix A and B) for the A,B polymer mixture is given by the well 
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known equations: 



and 



Xc " 2N A N B { ' 



& - Tm^M' (3 - 8) 

where xT^ * s the value of Flory-Huggins parameter at the critical point and is 
the critical concentration of A monomers (mf stands for mean-field). In order to 
include the fluctuations (i.e. dependence of \& on positions) we have to study the 
correlations functions in the system. Please note that the x parameter is rescaled 
by ksT. 

A. Vertex functions and critical point 

First of all, the critical point is marked by a very strong scattering, i.e. the 
divergence of the scattering intensity I(q) at q — > 0. This quantity is proportional 
to 

7(g) ~ (r 2 (q, -q))" 1 = 1 < *(q)*(-q) > . (3.9) 
Here the average is calculated according to the following formula: 

<...>=! J 

where Z and Z(*&) are defined in section II D. 

At the critical point, T 2 goes to zero for q — > 0. Additionally, 

< tf(q)tf(k)tf(-q-k) > 



r 3 (q,k, -q-k)=^ 



< *(q)*(-q) >< tf(k)tf(-k) >< tf(q + k)*(-q - k) > 

(3.10) 



vanishes at the critical point for q, k — > 0. The two equations, 

r 2 (o,o) = o, (3.ii) 

r 3 (0,0,0) = 0, (3.12) 
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locate the critical point in the (T, 0) plane. 

The method dealing with the calculation of the averages in the form given 
by Eq.(3.9) is described in detail in Ref.[6,7]. Here we use an expansion in the 
number of loops (Chapter 6, Ref.[7]), which in the case of a polymer mixture has 
great benefit of being an expansion in a small parameter proportional to 1/ \/N. 
In the limit of N — > oo, only the terms with no loops survive. The zeroth order 
approximation (no loops) corresponds to the random phase approximation 2 ) (see 
Appendix A and B, Eq.(B.ll) and the critical point location is given by Eqs.(3.7-8), 
thus it is the same result as the Flory-Huggins result. 

The equations for the vertex functions in the one-loop approximation are dis- 
cussed in Appendix C and Fig. 7 and Fig. 8 show the one loop diagrams which 
contribute to the vertex functions. Now we can proceed with the solution of 
Eqs.(3.11,12) in the one-loop approximation (Appendix C). We are interested only 
in the first corrections to the mean-field results. Thus, we assume the solution in 
the following form: 

Xc = x ™f(l + ai ), (3.13) 

and 

c = C'(l + a 2 ), (3.14) 

where ai and a 2 are small corrections of the order 1/ and X ? f and are given 
by Eqs.(3.7-8).In order to simplify the lengthy formulas, we also assume that the 
lengths of the monomers are equal, I a = Ib = I, but, of course, it is straighforward 
to perform the calculations without this assumption. Combining all these equations 
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we find Xc and 4> c in the following form: 



9v _ (v^+v^) 2 A , f J2\ , ( 4tt \ 



AT 



3/2 + jVj /2 \ //\ 3 / 27T \ 



x 



(a^)^-^(l) 3 ) 



(3.15) 



9c =-7==- — 7= 1 + 



f 36 \ (VnI-Vn^)VnI rr\ 
Wo/ 3 ; (v^+v^) 2 w 



, (N 3 A /2 - N 3 B /2 )VNX 3 

where A is the cutoff discussed in Appendix D. 
For A^4 = Nb = N, we have C = 1/2 and 



2 / / 12 
1+ - 



Xc = T7 1 + 



+ 



4n 



(3.17) 



3p l 3 C$VN 

where C = \[C\ + C 2 is a constant (see Appendix D). Assuming C = 1 and the 
volume fraction (the ratio of the volume occupied by monomers to the total volume) 
Pol 3 = 1 we find 



x ' = n{ 1 + {7n ]] - (3 ' 18) 



We can use this formula to estimate the shift of the critical temperature induced 
by fluctuations. Let T c m ^ = 400K for N = 10 3 . Then we find from Eq.(3.18) that 
T c = 319K. The shift, AT = T c m/ - T c is 81 degrees! Even for N = 10 4 we have 
still 8% correction to Xc/x™^- Of course, the constant C should be the fitting 
parameter in the analysis of experimental results and its value can be larger or 
smaller than one. Nonetheless, we think that it should be of the order of 1. 
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For Na > Nb, we find that C is larger than the mean-field value of the 
concentration of A- monomers, §™K For Na 3> Nb, it follows that 

For — > oo and A_b finite, the critical concentration approaches zero, while the 
critical temperature approaches T™* . 

Finally, let us discuss the order of the n-loop diagram in the limit of Na 3> Nb- 
In this limit (p c ~ V Nb / V Na and is very small. It enters the n-loop diagram as 
a multiplicative constant in the form l/(/) n+1 . Additionally, we have 1/q 2 contri- 
butions for n — 1 loops in the n-loop diagram. After grouping all these terms and 
using the formulas for the cutoff (Appendix D), we find that any n-loop diagram is 
proportional to I/^NaN^ 1 ^ 2 ), and it follows that the small parameter in the 
loop expansion is now 1 / \^N~b- Because of the 1 / ^Na factor we also find that in 
the limit of Na — > oo T c — > T™f even for finite Nb- This result is in accordance 
with de Gennes' predictions 1 -*, for the polymer in a solvent. As we can see, the 
loop expansion breaks down if either Na or Nb is of the order of unity, since then 
it is no longer true that the succesive n-loop terms are small in comparison with 
the (n-l)-loop terms. This analysis also prove that the loop expansion in polymer 
blends is indeed the expansion in small parameter, providing we chose the cutoff 
according to the prescription given in Appendix D. 

B. The Ginzburg criterion 

The critical temperature and concentration are the nonuniversal quantities de- 
pending on the details of the system. On the other hand, the values of the critical 
exponents are universal and do not depend on any microscopic details of the sys- 
tem. A set of critical exponents characterizes the whole class of critical points. For 
example the critical exponents for the Ising spin system, low-molecular mass bi- 
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nary mixtures and polymer binary mixtures are the same. All these critical points 
are said to belong to the Ising universality class 7 ''. As was mentioned in the In- 
troduction, the scattering intensity diverges on the approach to the critical point, 
according to the following formula: I(q — > 0) ~ (T/T c — 1)~ 7 . We know that for the 
Ising universality class 7 = 1.26, while in the mean-field theory we get 7 = 1. As 
we have seen in the previous section the fluctuations shift the critical temperature, 
however the most important effect of fluctuations is the change of the value of the 
critical exponents. In the experimental analysis 8- 13 ^ of I(q), in order to extract 
7, one plots log /(g) versus log(T/T c — 1). (see Fig. 9). It has been observed that 
away from the critical point the plot is a straight line with the slope of 45 degrees, 
indicating 7 = 1. As one approaches T c closer the straight line begins to curve. 
Finally, very close to the critical point it becomes a straight line again, but with 
the slope corresponding to 7 = 1.26. The change in the slope marks the breakdown 
of the mean-field theory. The criterion that tells us when this happens is called 
the Ginzburg criterion. The region around the critical point where 7 7^ 1 is called 
the Ginzburg region. Since I~ 1 (q) ~ r 2 ((/, —q), we can base our estimate of the 
Ginzburg region on the analysis of the two-body vertex function given by Eq.(3.9). 
Our strategy is to write it down as a sum of two terms: a mean-field one and a cor- 
rection to it. As we approach the critical point the corrections grow in comparison 
to the mean-field term. Finally they are larger than the mean-field term, indicating 
the breakdown of the mean-field theory. This happens inside the Ginzburg region. 

Close to the critical point and in the limit of q,p — > we can rewrite Eq.(3.9) 
and Eq.(3.10) in the following form: 

r 2 (o, 0) = A 1 (4>- 4>T f ) 2 + HxT f - x) + iW, x), (3.20) 
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and 

r 3 (0, 0, 0) = A 1 ($- + F 2 (0, x ), (3.21) 

where A x = 2(l/N A $™f) 3 + 1/N B (1 - (p™ f ) 3 ) and Fi and F 2 are the one-loop 
corrections to r 2 (0,0) and r 3 (0,0, 0), respectively. They have been discussed in 
detail in the Appendix C. Now we can express 4>™f and x™^ i* 1 terms of the critical 
concentration <p c and the Flory-Huggins parameter at the critical point Xc- Now 
Eq.(3.20) reads as follows: 

T 2 (0, 0) =A 1 $- 4> c f + 2( X c - x) + (*i(0, x) - Fitfc, Xc)) - 2F 2 (^ C , X c){4> - 4>c)- 

(3.22) 

This equation is the basis for the quantitative formulation of the Ginzburg criterion, 
namely, the mean-field theory is valid providing the first two terms in Eq.(3.22) 
(giving the mean field result for 7) are much larger than the last two terms , which 
are the corrections. Thus the mean-field theory breaks down when 

Ai{4> - 4> c ? + 2( X c ~X)< ~(Fi(4>, x) ~ F^c, Xc)) + 2F 2 (^ C , X c){4> - 4>c). (3.23) 

Now, we can give the explicit formula for the Ginzburg criterion in the special case 
of Na ~ Nb and 4> ~ 4> c . Neglecting the dependence of I4 on k, i.e. assuming 
the approximation T^^k, — k, 0, 0) « T^^O, 0, 0, 0) we find the Ginzburg criterion 
in following form: 

2(Xc _ x) + Al{$ _ h? = £f^±S^r^ . (3 . 24 ) 

p 2 {R A N A (j) c + R B N B (1- 0c) 1 ) 3 

Here C g is a constant, Rf = iV;/ 2 /6 (i = A,B). For = C , Eq.(3.24) reduces to 
Eq.(3) of Ref[10] providing we change v m to l/po, 4>™f to C , xT^ to Xc and C 
to Cg/2. On Fig. 10 we show the plot of the scattering intensity versus the inverse 
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of the temperature 10 ^. Please note that the mean field critical temperature is not 
correctly identified on this plot 15 ) (see also 74 )), since we find that 



I(q - 0) ~ ( Xc - x) 



-l 



(3.25) 



where \ c is not the mean-field Flory Huggins critical parameter. In fact 
the Ginzburg region is of the order of 1/N whereas the difference between the 
critical temperature and the its mean field value goes as l/y/N. Therefore it is 
not correct to identify the mean field critical temperature with the point where the 
extrapolated straight line representing I(q — > 0) far from the critical point crosses 
the temperature axis. 

Finally, we consider the limiting case of Na ^> Nb- In the first approximation 
and for = C , one finds that the mean-field theory breaks down for 



Thus, for Na — > 00, one would expect that the mean-field theory is valid in the 



whole region around the critical point, for any value of N B (Note that l/y/N a 

term here comes from the l/y/N a term which apppears in the loop expansion). 
This result is misleading for the following reasons: In our analysis we have not 
specified the temperature below which we can expect any influence of the critical 
point on the behavior of the polymer mixture. Such a temperature is called the theta 
temperature and is denoted, O. Above this temperature mixing is strongly favored 
and the influence of the critical point on the scattering intensity is negligible. Only 
below can we expect any sign of critical behavior (either mean- field or Ising-like). 
This means that if the Ginzburg region estimated above is larger than the region 
between O and T c we will never see the mean- field behavior, characterized by 7 = 1, 
but only the Ising-like behavior, characterized by 7 = 1.26. As shown by Joanny 5 ) 



(Xc - x) 



1 



(3.26) 



Xc 




T-T c 1 



(3.27) 



®-T c y/N B ' 
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marks the breakdown of the mean-field theory in the aforementioned sense. Finally 
we see that in the limit of N B — > 1 the mean-field theory breaks down completely 
even for infinite Na- As mentioned in the previous section also the loop expansion 
breaks down, but peculiarly the critical temperature approaches its mean-field value. 

C. Single chain statistics 

The single chain statistics is determined in the first place by the correlation 
function for the single chain. In the A,B homopolymer mixture we chose one chain 
(e.g. A) and tag it by coupling it to the fictitious field Ua- The correlation function 
is given by 17 ^: 



w qi ,q 2 ) =< &x*m>*) >= z[u JuZ)Iu aM - (3 ' 28) 

where the partition function (section II) Z as a functional of the external field Ua 
is given by 

Z[U A ] = J D^ A) J ^4 nB) exp(-^[A nA) ,4 nB) ]) 

J DJ A J DJ B <*p(i J J^^JA + i J J^A b) Jb (3.29) 

+F A nA) [J A ] + f£ ] [J a + iU A ] - [ J a] + F B nB) [ J B ] j . 

Here F^ 1 ^ is defined in section II and Appendix A. Now using this equation at 
Ua = we find the following approximate equation for the single chain structure 
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factor in the incompressible blend: 



VS AA (d, -q) =< ^(q^C-q) >0 + I 

z Po J 

- < ^(q^C-q) >o< ^(k^-k) > ^ 

r^(k)r^(k) 



dk 



(2tt) 3 



(3.30) 



r 2 (k) 



It includes the ideal term and the first nonvanishing correction to it. Here 1?2 is 
defined in Appendix C and in Appendix A. As defined in the Appendix A the 
averages with the subscript are taken over the configurations of noninteracting 
system (ideal averages) whereas those without the subscript are taken over the 
configurations of the interacting system. 

The equation for the single chain correlation function in Ref[20] is similar in 
structure to Eq(3.30), only instead of having r^ ) (k)r^ ) (k)r" 1 (k) one has there 
< Jyi(k) Jyi(-k) >. The latter is exactly equal to %(1 + 2p xT2 1 (k))/2. Note that 
if we use the linear approximate relation between J a and ^a i-e. J a = r^^A 
we would get the same result from both approaches. This equation for the single 
chain correlation function should be supplemented by the self- consistent equations 
for the collective structure factor, T 2 (k) (Appendix C). The equation for the single 
chain correlation function and the collective vertex functions are also coupled via 
the radius of gyration (see Appendix D) given by 75 - ) : 



It is obtained from Eq(3.30) by differentiating both sides twice with respect to q 
and taking the limit of q — > 0. We find 



2^/3 = 2NlR 2 /3 - -L y ^f(k)4 A \k)4 A \k)T^(k) (3.32) 




N A N A 



(3.31) 
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where Rq = Na1 2 /6 is the radius of gyration for the ideal chain (Eq(3.3) with the 
ideal average) and /(k) is explicitely given in Appendix E. Since f(k) is always 
positive we find the first intuitively obvious result that chains in a blend shrink 
below the Gaussian chain limit. 

The role of the cutoff (Appendix D) is also clear here. If not for the \fN in the 
upper wavevector cutoff the correction to the Gaussian radius of gyration would 
be nonzero even in the limit of N — > oo and thus RPA would not be a correct 
description in this limit. 

IV. Linear flexible diblock copolymer system 

An AB diblock copolymer molecule consists of two blocks (A and B) joint 
together. The blocks are incompatible, but contrary to the case of AB homopolymer 
blend cannot separate on the macroscopic scale due to the chemical bond joining 
them. Nonetheless as we lower the temperature the blocks separate on the scale of 
the radius of gyration forming the periodic structures - the phase transition between 
the disordered phase and ordered is called, in the polymer community jargon, the 
microphase separation. In the case of the diblock copolymer the concentration of A 
monomers is set by the architecture of the chain, the average fraction of A monomers 
is given by 

t NA (A<i\ 

f= W , (4-1) 

where N = Na + Nb is the total number of monomers in the diblock copolymer 
molecule. The order parameter field, 

*(r) = 0j ) (r)-/, (4.2) 



is defined as in section II and n is the number of copolymers. Incompressibility 
constraint is implied. Contrary to the case of homopolymer blends the equilibrium 

28 



distribution of A monomers at low temperatures, given by is not uniform in 
space. 

A. Scattering intensity in the random phase approximation 

Leibler 22 ) applied random phase approximation (RPA) scheme (described in 
Appendix A and B, see also Ref[76,20]) to the system and found the following form 
of the two body correlation function (Eq(3.9)): 

Wi(q) 

I(q) = W 2 (q)-2p oX W 1 (q) (43) 

where W\(q) = S11S22 — Sf 2 and W 2 (q) = Su + S22 + 2Si 2 depend only on the 
architecture of the diblock copolymer (see section II) ; W\ is the determinat and W 2 
the sum of all the elements of the matrix S a p, where 

S11 = ly < 4%\q)<lft\-q) >o= -<7i (/,*), (4.4) 



v -<<P ( B } m ( B\-Q) >0=^ 

v Po 



S22 = ^< ^\q)^\-q) > = -gi(l - f,x) (4.5) 



and 



S12 = S 2 i = ^-< ^\q)^\-q) >o= ^~ (9i(l,x)- gi (f,x)- gi (l - f,x)) 

(4.6) 

Here x = q 2 Nl 2 /6 : V is the volume and 

9l (f, x) = 2(fx + exp(-fx) - l)/x 2 . (4.7) 

Please note (compare with Appendix A) that since the blocks A and B are joined 
the cross corelation function S\ 2 is non zero here. The behavior of I(q) is shown 
in Fig.ll. As we increase x (lower the temperature) the peak of /(g) grows and 
finally diverges at the well defined (and fixed) q* vector. Interestingly the position 
of the peak in this approach is fixed independent of the temperature. For f=l/2 
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I(q) diverges at xN = 10.55, while for / = 0.25 at xiV = 17.6. The divergence of 
the scattering intensity indicates the onset of ordering at the typical length scale 
given by the inverse of the scattering vector q*. 

The order parameter depends on position in the spatially ordered phase. The 
free energy in the random phase approximation has the following form in this case: 

= E / • • • ^ys r ^ • • - q»)*(q0 • • • *(q»)> ( 4 -«) 

where the vertex functions are given in Ref[20,76,22] and the general method 
of the calculation has been described in Appendix A and B. The equilibrium con- 
figuration of the system corresponds to the minimum of O i.e. to the solution of 
the functional equation: 

In practice this complicated equation is simplified if we assume a given form (in 
accordance with the symmetry of the phase) of the spatial dependence of \&(r) in 
the ordered phase. For example in the lamellar phase: 

*(r) = Acos(q t z). (4.10) 

Now inserting (4.10) into Eq(4.8), we find the free energy in the lamellar phase as 
a function of two variables, the amplitude, A and the wavevectors qi. The period 
of the lamellar phase is equal to 2n/qi. Now the minimum of the free energy in the 
lamellar phase corresponds to the solution of two algebraic equations: 

dQ(A, qi ) 



OA 

and 



= (4.11) 



= o. (4.12) 

dqi 
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The phase transition between the disordered phase and the lamellar phase takes 
place at the point where both energies are equal. In the random phase approxi- 
mation the transition is second order and takes place for qi = q* and at the point 
where the scattering intensity diverges. Apart from the lamellar phase one studies 
in the same way 22 ), other phases such as hexagonal phase, bcc phase, obdd phase 
or others. For example the simplest approximation of \1/ in the hexagonal phase has 
the following form: 

3 

vl/(r) = A^cos(q n r ± ), (4.13) 

n=l 

where q\ = (1, l/V^)^, q2 = (— 1, l/y/3)qh and qs = (0, 2/y/3)qh are the vectors 
spanning the first shell (similarly in the lamellar phase Eq(4.10)) in the reciprocal 
space for the regular hexagonal lattice. The form of \& given by Eq(4.10) for the 
lamellar phase and Eq(4.13) for the hexagonal phase is a good approximation close 
to the weakly first order or second order phase transitions, it breaks down at low 
temperatures, when the ordered domains are sharp 76 ). 

In the random phase approximation the hexagonal phase and bcc phases are 
also stable in some range of parameters. The transitions between these phases and 
disordered phase are first order. The phase diagram of the diblock copolymer system 
is shown in Fig. 12. In the weak segregation regime the random phase approximation 
is exact in the limit of N — > oo. 

B. Fluctuations induced first order phase transition 

The transition to the lamellar phase from the disordered phase is second order 
in the random phase approximation. It can be easily seen if we use Eq(4.10), apply 
it to Eq(4.8) and expand the free energy in the amplitude A. We find the following 
form of O (keeping the first three terms only): 

O = a A 2 + bA 4 + cA 6 (4.14) 
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Here b and c are two positive constants independent of the temperature and a 
changes sign at the point when the scattering intensity I(q) (Eq(4.3)) diverges. At 
this point we have the second order phase transition. Now if we go beyond the ran- 
dom phase approximation and include the one-loop corrections the second term in 
Eq(4.14) becomes negative before a does, but c remains positive and consequently 
the transition becomes first order 20 ' 23 ' 77 ). The free energy in the one loop approx- 
imation has the same general form given by Eq(4.8), but the coefficients in the 
expansion are given by the vertex functions T n calculated in the one loop approx- 
imation. The equation for the two and three body vertex function in the one-loop 
approximation is given in Appendix D. These equations are general and valid for 
polymer blend system as well as for the disordered phase of diblock copolymers. In 
the ordered lamellar phase the equation for vertex function are slightly different, 
since the field ^ given by Eq(4.10) is non-uniform. In the case of/ = l/2r 3 = 
due to the symmetry 22 ) and the one-loop equation for the lamellar phase has the 
form given in Appendix F 20 ). The phase diagram in the one- loop aproximation is 
shown in Fig. 13. As we can see the transition to the lamellar phase occurs not only 
for the symetric case / = 1/2 (as we have found in the random phase approxima- 
tion), but also in the neighbourhood of this point. The transition to the lamellar 
phase is first order as mentioned previously. In the limit of N — > oo the one-loop 
corrections vanish and the results reduce to the random phase approximation. 

In the random phase approximation the peak in the scattering intensity is fixed 
and scales with the polimerization index as \/~N. However from the experiments 78 ) 
it follows that the peak in the structure factor exibits a much stronger dependence 
on N than predicted by RPA. In Fig. 14 the comparison of the results of the one-loop 
approximation 20 ) and experiment are shown. The logarithm of the peak position 
is plotted against the logarithm of N. The agreement between the theory and ex- 
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periment is remarkable, especially that no adjustable parameters has been used. 
The vertical straight line is the predicted N for which the microphase separation 
occurs (the temperature is fixed here at T=296 K). Close to this transition the peak 
position q* is approximately represented by N~ - 8 . 

Since the peak position deviates strongly from the V~N behavior, it indicates 
that the conformation of a single polymer molecule should also differ from the Gaus- 
sian behavior. Indeed, both in the theory 20 ) and computer simulations 21 ), strong 
deviation from the Gaussian behavior has been found. The equations for the radius 
of gyration for the polymer blend are discussed in section III, the very similar equa- 
tions are found for the diblocks 20 ). The solution of these equations indicate that 
the chains strongly stretch as the temperature is lowered towards the microphase 
separation transition temperature. The configuration of the chains indicates that 
stretching occurs at the point where the two blocks are joined. However the blocks 
themselves shrink below the gaussian limit. Thus the configuration of the chain re- 
sembles the dumbell, consisting of two coils more densely packed than the Gaussian 
coils, and stretched above the Gaussian limit. The fact that the coils shrink below 
the Gaussian limit has the same origin as the shrinking of the polymer chains in the 
blends close to the critical point (section III C). Interestingly the deviations from 
the Gaussian behavior occur far from the microphase separation temperature, deep 
in the disordered phase. 

C. Gyroid phase and self consistent field theory 

The diblock copolymer system forms a variety of phases. So far we have only 
mentioned the simplest structures studied in the theory: the lamellar, hexagonal 
and body center cubic structures. Apart from those simple morphologies the diblock 
copolymer system exhibit new type of structures: the bicontinous phases where 
the junction between the blocks are located at the periodic continuous internal 
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interfaces. The double diamond and gyroid structures are shown in Fig.l and 2. 
Both were observed in the polystyrene-polisoprene system. The first one in the 
strong segregation regime 30-32 -' (which is not discussed in this paper), while the 
second one in the weak segregation regime 79 ^, where the Landau theory should 
be valid. This phase has not been studied in the random phase approximation, 
but in the mean- field theory called the self consistent field theory 24 ). The latter 
theory applied to the system of diblock copolymers predicts the stable gyroid phase 
between the lamellar and hexagonal phase as shown on Fig. 15. In this theory 
the problem of computing the partition function for n chains is reduced to the 
single chain in the external field. The difference between this approach and the 
random phase approximation starts at the level of the partition function for the 
system of noninteracting chains in the external fields J a, Jb (Appendix A Eqs(A.7- 
8)). In the calculations this function is computed exactly in the self consistent field 
theory, while in the random phase approximation it is expanded in Both theories 
neglect fluctuations. The gyroid phase ha not been studied in the random phase 
approximation. 

The gyroid structure is stable for %iV between 14 and 20, which is in agreement 
with experiment. Please note that for %iV larger than 20 the double diamond 
structure is stable (see Fig. 3). There is no direct transition between the disordered 
and gyroid structures. Instead at / = 0.452 and xN = H-14 there is a triple point, 
where the lamellar, hexagonal and gyroid phase coexist. Below %iV = 11.4 the 
hexagonal phase transforms directly into the lamellar phase. In this approach the 
double diamond (obdd) structure is only metastable, but as we mentioned previously 
it is understable since the theory is applied to the weak segregation regime, while 
the obdd phase is expected at large xN m the strong segregation regime. 

The same method has been applied to the diblock copolymer system with 
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conformational asymmetry, i.e. to the system where the length of the monomers 
/ is different for different types of monomers. The phase diagram in the case of 
Ia/Ib = a/TO is show in Fig. 16. By comparing the diagram in Fig. 15 (Ia/Ib = 1) 
and Fig. 16 we note that conformational asymmetry does not change the topology 
of the phase diagram 25 ) . 

V. Other flexible block copolymer systems 

The variety of different block copolymer architectures is infinite. Therefore, 
here we shall chose only a few, to prove that the methods discussed so far can also 
be applied to block copolymers irrespective of the specific architecture. Here we 
discuss, in the random phase approximation, the triblock copolymers, diblock ring 
copolymers and random copolymers. 

A. Ring copolymers 

The AB diblock ring copolymer consists of two blocks joined in two points. The 
ring architeture can be easily included in the general scheme presented in section 
II. The same method as applied to polymer blends or linear diblock copolymers and 
described in Appendix A, B and C can be also applied to diblock ring copolymers. 
Here we consider the scattering intensity and the instability (same as in section IV 
A) in the system. 

The scattering intensity I(q) for ring copolymers, calculated in RPA, has ex- 
actly the same formal structure as given for linear diblock copolymers in section IV 
A (Eq.(4.3)). Only the matrix S a p is different 48 ) , i.e. 

Su = >L f ds [ f cfa'e<-*l— 'IKM— 'D, (5.1) 
Po Jo Jo 

S 22 = ^ f'ds f 1 cfa/e^l— 'D(i-I— 'D, (5.2) 

Po Jf Jf 

S 12 = - [ f ds t cfa'e'-*"— 'IK 1 "!— '!>, (5.3) 
Po Jo Jf 
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x is defined after Eq(4.7). The same formulas are obtained for the linear diblock 
copolymers but without the term 1 — \s — s'\ in the integrand. Due to this term the 
integrals cannot be calculated analytically for ring copolymers. 

Let us first compare the scattering intensity of the linear and ring diblock 
copolymers in the case of x = (high temperature) and / = 1/2 (symmetric case). 
In Fig. 17 the scattering intensities are shown. The peak in the ring copolymer is 
at larger q vector (which could be expected) and is less intense than the peak for 
the liner diblocks. In Fig. 18 the instability line is shown (at which the scattering 
intensity diverges) . As we see the rings require a larger amount of repulsion between 
A and B monomers to segregate. The latter result is also understandable since the 
block in the ring are more constrained than in the diblock. To see how severe are 
these constraints it is instructive to compare the instability in the ring copolymer 
to that in the (AB) n copolymer 80 ) where there are n AB diblocks joined together in 
the linear chain. The size of the AB diblock is N and each blocks has the size N/2 
(symmetric case) . It appears that in the limit of n — > oo the critical Flory Huggins 
parameter approaches a limiting value %A — > 15, whereas for the symmetric diblock 
rings (see Fig. 18) of size N the instability occurs at %A = 17-8, thus the transition 
occurs at even lower temperature. In the diblocks the same instability occurs at 
xN = 10.5, thus the transition temperature to the lamellar phase is approximately 
40 % lower in the ring system than in the linear diblock copolymer system. 

The result for the microphase separation obtained in the random phase ap- 
proximation has been confirmed in the computer simulations 47 -' . This result is in 
fact surprising, since in the random phase approximation it is assumed that chains 
adopt the Gaussian configurations, whereas in the ring system the constraints are 
such that the statistics is not Gaussian even in the limit of N — > oo. The scaling 
exponent for the radius of gyration is v = 0.45 (R g ~ N u ) as has been found in 
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computer simulations 47 ' 81 ' 82 and theory 46b ). The comparison of the radius of gyra- 
tion for the linear diblock chain and ring diblock chain is shown in Fig. 19. Here T c 
is the transiton temperature to the lamellar phase and roughly concides with the 
instability temperature shown in Fig. 18. 

We have learned that the linear diblock copolymer system stretch as the tem- 
perature is lowered 20 ' 21 ^. The same applies to ring diblock copolymers 47 ) . As the 
temperature is lowered the strong increase (20-25%) in the distance between the 
center of masses of two blocks is observed (Fig. 20). At the same time the blocks 
themselves shrink (Fig. 21) similarly as in the case of linear diblock copolymer sys- 
tem discussed in section IV. The effect of the stretching is more pronounced in the 
rings (20-25% above the radius of gyration at the infinite temperature) than in the 
linear diblocks (15% above the high temperature radius of gyration). 

B. Triblock copolymers 

The ABA triblock copolymer system has been studied by several authors 83-86 ) 
The phase diagram for this system is similar to the one observed in diblock copoly- 
merg 84,85) (pig 22). This observation has been supported by experiments 83 ' 87-89 ). 
This is not surprising since cutting the triblock in the middle produces two diblock 
coplymers and the entropy gained from doing this should not be very large, so the 
free energy of triblocks in a particular phase should be roughy the same as that of 
twice as many diblocks of half the length in their analogous phase. 

The scattering intensity for the triblock copolymer system has the same gen- 
eral form as given by Eq(4.3), only the matrix S a p is different 83 -* (gi is defined in 
Eq(4.7))i.e. 

Sn =- (<7i(/i, N) + 9l (h, N) + 9l (f 3 , N) 

P° V , (5.4) 

+<7!(1, N) - 9l (l - / 3 , N) - 9l (l - f u A)) 
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S22 = —gi(h,N), (5.5) 
Po 

S12 ~ (gi(l - fi, N) + 9l (l- f 3 , N) + 9l (l, N) 

2 P° V , (5.6) 

-gi(fc,N)-g 1 (f 1 ,N)-2g 1 (f 2 ,N)} 

Here /i, f 2 , fs are the monomer fractions in each block, 1, 3 are A blocks. In the case 
of fs = (or fi = 0) the formula reduces to the diblock copolymer formula given 
in Eqs(4.4-4.6). 

We can also expect that as we lower the temperature the chains stretch, but 
the individual blocks shrink. 

There is one notable difference between the system of ABA triblocks and AB 
diblocks. In the former case there are bridge and loop configurations (Fig. 4), de- 
pending on the location, in lamellas, of the different A blocks belonging to the same 
triblock. The presences of bridges should affect the mechanical properties of the 
triblock system 86 ) or in general the multi-block copolymer system 44 -*. It has been 
found that about 40-50 % of copolymers form bridges in the triblock copolymer 
lamellar phase. 

C. Random copolymers 

In the AB random copolymer system 40 ) the architecture of the chain is specified 
not only by the fraction of A monomers in the chain, /, but also by the conditional 
probabilities that the A monomer in the chain is followed by the A, paa, or B 
monomer, pab and similarly for B monomers. One naturaly has paa + Pab = 1 
and same for B. The function / and A = paa + Pbb — 1 are sufficient to describe 
the disorder quenched into an AB random copolymer. Now apart from the densities 
already introduced in section II we have to include the occupation number 9i at the 
i-th point along the chains. As in all the systems with quenched disorder one has to 
calculate all the relevant thermodynamic quantities for fixed distribution of 9 and 
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then average them with respect to 6>, which describes the disorder. The practical 
way to do it is by the replica method 40-43 ' 90 ^. 

Let us study the macrophase separation in the melt, which is similar to the 
separation in A, B homopolymer blend. The condition for the instability in the 
melt in the random phase approximation is: 

4 /Yl-fW-/= 2N(1-X) 2 

J{ J)Xc N(-X) 2 + 2X(N(1-X) + X N -1) K ' 

For A — > 1 this condition reduces to the same equation as we would get for the 
mixture of A homopolymers and B homopolymer each of length Nl i.e. 

4/(l-/)xr / = |- (5-8) 

In the limit of A — > we would have a system where the copolymers have the very 
short sequences of the A and B monomers. Thus the condition (5.7) reduces to 

4/(1 - f)xT f = 2 (5.9) 

as for the system of disconnected monomers. As we see in this limit the interesting 
scale of the phenomena is related to the short sequence of the A or B monomers 
and not to the size of the polymer molecule. Consequently our description is no 
longer valid; in particular the structure and phase transitions in the system depend 
on the details of the intermolecular potential, which in the present treatment have 
been neglected. 

VI. Mixtures of homopolymers and diblock copolymers 

So far we have considered the polymer melts or blends which were described by 
a single order parameter. In order to illustrate the application of the method pre- 
sented here to more complicated polymer blends we consider a ternary mixture of A, 
B homopolymers and AB diblock copolymers 34-37 -*. This system is in some respects 
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similar to the ternary mixture of oil water and amphiphiles 38 ' 91 ''. The copolymer 
simply acts as a surfactant accumulating at the interface of A-homopolymer rich 
phase and B-homopolymer rich phase, reducing in this way the number of unfa- 
vorable contacts between A and B monomers 36 ' 92 ^. This analogy will be exploited 
here. 

A. Structure factor matrix in the RPA 

To simplify the calculations we assume that the length of all the polymers is N 
and the diblocks are symmetric. It is convenient to introduce a set of three order 
parameters. The first order parameter, rj, measures the deviation from its average 
value of the concentration difference of A and B homopolymer. The second e, 
measures the local deviation from its average value of the copolymer concentration. 
The third, measures the local deviation from its average value of the total fraction 
of A-monomers. All order parameters vanish in a spatially uniform, one phase 
region. 

Now the scattering intensity 

%) = (r^ 0) ( ? ))" 1 (6.1) 

becomes a matrix and so does the vertex function . For the matrix 



r (0) 
1 2 



J- 7777 




r?7* 


r 















(6.2) 



calculated in the random phase approximation (Appendix A,B) we obtain 
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'> = -( 75 I (C5) 

&AA — JAB 

2 

r** = 2x, (6.7) 



where 



5 A = ^-^ Jv ^( 1 ' a; )' ( 6 - 8 ) 
fif B = ^l^^ 1 (l,x), (6.9) 

Saa = —N 9i (1/2,x), (6.10) 
Po 

•^ab = ^N( 9l (l, x) - 2 9l (l/2, x)), (6.11) 

and </» is the average concentration of copolymer. For simplicity we have assumed 
that the concentrations of homopolymers A and homopolymers B are equal. Func- 
tion <7i is given by Eq(4.7). 

The two most useful scattering intensities, are obtained by taking the inverse 
of the vertex function matrix Eq(6.1), and are given below: 

(Sa + Saa — Sab) 



2[1 - Pox(Sa + S A a - S A b)] 
in agreement with Broseta and Fredrickson 34 ) and 



(6.12) 



T 2S a [1-x(Saa-Sab)] 

m ~ [1-Pox(Sa + S A a-S ab )Y [bA6) 
in agreement with Leibler 92 ). 

For completeness, we compute also the copolymer copolymer structure function 
at equal homopolymer concentrations, 

ee ^ Sa Saa + Sab^ 

= N gi (l,x)<f>(l-<f>)/po. (6.14) 
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This result, in agreement with Leibler 92 ), is independent of temperature, and al- 
ways has a peak only at zero wavevector only. In this respect it is similar to 
calculated amphiphile-amphiphile structure functions in systems of oil, water and 
surfactant 38 ' 94 ). 

B. Stability limits, disorder, Lifshitz and equimaxima lines 

Much of the phase diagram can be obtained from either of the above structure 
functions. The limit of stability of the disordered phase occurs at the highest 
temperature at which a divergence in the above structure functions appears. From 
Eqs. (6.12) or (6.13), we see that this happens at the smallest value of x- Here we 
treat x as a function of x so that we can minimize it with respect to x. From the 
denominator of the structure functions we find 

' =hs A (x) + S AA (x)-S AB (x)), (6.15) 



p oX (x)N N 
or, using the explicit formulas, 

1 e" x (l -2(f)) +4(/)e- x / 2 +x - 1 -20 



x(x)N x 



2 



(6.16) 



where again, x = q 2 R 2 = q 2 N 2 l 2 /6. Minimizing x(x), we find that for < 2/3, the 
minimum occurs at x = 0, and that its value is 

XN =T ^, 0<2/3. (6.17) 

Along this line, the disordered phase makes a continuous transition to uniform 
(q = 0) phases. This is just the macroscopic phase separation into A-rich and B- 
rich phases. For a copolymer density between 2/3 and 1, the instability occurs at a 
non-zero value of x and the continuous transition is to a phase characterized by the 
corresponding non-zero wavevector (x) 1 / 2 /R; this is the microphase separation to 
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a lamellar phase. The density of copolymer at which a transition occurs to a phase 
characterized by x 7^ is 

-2 + x + (x + 2)e- x 
4>{x) = — T -, — N , — x 7^7, < x < 3.785 6.18 

and the value of %iV at which it occurs is given by Eqs. (6.16) with x in the range 
indicated in Eq. (6.18). The transition lines are shown in Fig. 23. There is a region 
of three phase coexistence between the A-rich, B-rich and lamellar phases which is 
shown schematically in the figure. 

Much more information about the disordered phase can be extracted from the 
structure functions, of course. To obtain some information on how the copolymer 
tends to order the homopolymer and also to order itself, it is useful to introduce the 
disorder line and the Lifshitz lines 35 ' 38 ' 39 ). The disorder line is the locus of points in 
the phase diagram at which a general correlation function, G(r) : no longer decays 
monotonically with distance at large r, but acquires a component which, at large 
distances, behaves like an exponentially decaying oscillatory function; 

G( r ) -> ex P(~ r /£) S in(27rr/A + 5). (6.19) 
r 

This oscillatory behavior reflects the tendency of the copolymer to order the A and 
B monomers in space. Thus, on the copolymer-rich side of the disorder line, two 
lengths are needed to describe the disordered phase; the usual correlation length £, 
and a wavelength A. As the amount of copolymer in the system is increased, one 
expects the wavelength to decrease as order can be enforced over shorter distances; 
conversely, as the copolymer density decreases, A increases and diverges at the 
disorder line. We locate the disorder line as follows. The Fourier transform of a 
general correlation function can be written as a linear combination of the Iij(q), 
and the poles of these functions determine the behavior of the correlation function. 
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If these poles are located at a value of q which is solely imaginary, the correlation 
function decays monotonically with distance. Thus the disorder line is the locus 
at which these poles first acquire a non-zero real part. As 7y = (I2 j" 1 , the 
location of these poles is also the locus at which the determinant jr^ ^)! vanishes 
for imaginary q. Because this determinant is common to all iy, all correlation 
functions have the same disorder line in general. From the structure functions 
Iw or I vv we see that the disorder line is found by locating the poles of either 
correlation function, just as the phase boundary was located. However to locate 
the phase boundary, we sought the poles occurring for real q (corresponding to 
x > 0), whereas to locate the disorder line we seek the poles occurring for q just 
as it becomes purely imaginary (corresponding to x < 0). We do this by setting 
x = —y + id, y > in Eq. (6.16), and equating real and imaginary parts. Then 
taking the limit of 5 to zero, we obtain two equations which can be solved for <f>(y) 
and x(y)N on the disorder line. They are 

^ = 2[ev(2-y)-ev/H4-y) + 2Y ^ 



and 

1 e 3y / 2 -e y (y + l) + e y / 2 (y-l) + l 



(6.21) 



X (y)N y[ e y/2^_ y) _ 2 - e y(2-y)] ' 
where y ranges over all positive real numbers. Two limits readily emerge from 
these equations. They are that <p approaches 2/3 as xiV approaches 6, and that <p 
approaches 1/2 as %iV vanishes. The disorder line is shown in the phase diagram 
of Fig.24 as the line LD. 

The disorder line shows where oscillatory behavior first appears in the corre- 
lation functions of the disordered system. The locus at which this behavior first 
dominates a particular correlation function is given by one of two kinds of line. The 
first is the Lifshitz line; at this line, the peak in the associated structure function 
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moves continuously from zero to a non-zero wavevector. Its location is determined 



In general, each structure function has its own, distinct, Lifshitz line, and the 
difference in these lines provides information on the different effects the copolymer 
has on the various components of the system. The Lifshitz line associated with 
the structure function of all A monomers, is readily found to be given by = 2/3, 
independent of temperature. This line is shown as LF in Fig. 24. The fact that the 
disorder line and the Lifshitz line of Iy-q, are close indicates that once the tendency 
of the copolymers to order the A monomers appears in the correlation functions, 
relatively few additional copolymers are needed to make this tendency dominant 
in the total A, total A correlation function. It is not clear, however, whether the 
copolymer affects the A monomers in the copolymers or in the homopolymers more 
strongly, or whether it affects them equally. This can be determined by examining 
Irjrj, which depends only on the homopolymer concentration, and finding the locus 
at which this function is dominated by components of non-zero wavenumber. In 
the vicinity of the point xiV = 6, (p = 2/3, this dominance is seen as the peak at 
zero wavevector moves to nonzero q, a movement which occurs at the Lifshitz line 
of this structure function; 



However, another mechanism takes over beyond <fi « 0.79. One finds instead that 
a second peak at non-zero q develops in I vv , and becomes larger than the peak at 
q = 0. This occurs before the Lifshitz line condition, Eq. (6.22), is satisfied. This 
new line, which we refer to as an equimaxima line, is defined as the locus of points 
at which the two peaks, one at zero, the other at non-zero wavenumber, are of equal 



by 



«9%-(q) i 

dq 2 '<?=o 



= 0. 



(6.22) 




(6.23) 
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height. At this line, the wavevector characterizing which is the largest peak in the 
structure function moves disco ntinuously from zero. The line is determined by the 
simultaneous solution of 

dl m {x,(j), x) I 



dx 

and 



0, (6.24), 



I V r,{x*, 0, x) = W°> 0> X)- (6-25) 
The equimaxima line (EG on Fig. 24) and Lifshitz line (LE) join smoothly as x* — > 
which occurs at (p = 121/153 0.79 (point E). We clearly see from Fig. 24 that even 
though oscillatory components enter all correlation functions at the disorder line 
and dominate the correlation function of all A monomers with only a rather small 
additional increase of copolymer concentration, they do not dominate the correlation 
function of those monomers only in the homopolymers until a large increase of 
copolymer concentration, so large that the system is about to become unstable to 
the lamellar phase. We conclude, therefore, that the copolymer is very efficient in 
self organization, but very inefficient in the organization of homopolymers. 
C. The Landau-Ginzburg free energy and the Lifshitz point 
To simplify the calculations we consider the case of equal homopolymer concen- 
trations and integarte out the concentration of copolymers e. Therefore we are left 
with one single order parameter r\. We find a Landau free energy in the following 
form: 



^k B Tp J dr ^h \ 



Q[rj] = -k B T Po I dr[ h \V 2 V (r)\ 2 + go\V V (r)\ 2 + f(rj) ] (6.26) 



where the bulk free energy is: 

f{rj) = arj 2 + bif + of. (6.27) 

Here go and ho are constants indepent on r), but in general they do depend 36 - 1 on 
r\ ih[r\ = 0) = ho and g(r] = 0) = go). The parameters a, 6, c and ho, go are given 
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in the Appendix G. An ordinary critical point occurs when a = 0, while the other 
parameters are positive, and an ordinary tricritical point occurs when a = and 
6 = 0, with the other parameters positive. Because go is positive, the ordered 
states which occur just below these transitions are spatially uniform, characterized 
by a vanishing wavevector q. At a Lifshitz point, the stability of these uniform 
states is lost because g = there, so that these points mark the last appearance of 
uniform ordered states and the first appearance of non-uniform ordered states (e.g. 
a lamellar one) characterized by a non-zero wavevector. At a Lifshitz critical point, 
a = 9o = 0; similarly at a Lifshitz tricritical point, a = b = go = 0, with all other 
parameters positive. One sees immediately how unusual this point is expected to be 
as three coefficients must manage to vanish simultaneously. As you can see from the 
results given in Appendix G the point L on Fig. 24 is the Lifshitz tricritical point. 

We first determine the behavior of the uniform solution, r]o, as the Lifshitz 
tricritical point located at T c is approached (go = b = 0, a ~ t = (T — T c )/T c ) by 
variation of the free energy of Eq(l.l) and find ijo ~ with f3 = 1/4. Inserting this 
into Eq. (6.26), we find F ~ t 2 ~ a with a = 1/2. These exponents also characterize 
the ordinary tricritical point in mean field theory. Thus the Lifshitz nature of 
the tricritical point does not affect these exponents within mean field theory. To 
see the difference between the Lifshitz and the ordinary tricritical point, we look 
at the behavior of deviations of the order parameter from uniformity. We find 
that r^ 0) (q -> 0) ~ V with 7 = 1, again independent of the Lifshitz phenomena. 
However, if we rescale the wavevector q by the correlation length £, we find that 
£ ~ t~ u with v = 1/4 in contrast to ordinary critical behavior in which, with 
go nonzero, v would be 1/2. The Lifshitz behavior also is manifest in the value 
of the upper critical dimension, d*, the smallest dimension at which mean field 
behavior is correct. For ordinary tricritical behavior, d* = 3, so that fluctuations 
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are negligible. (They contribute logarithmic corrections to mean field results.) For 
Lifshitz tricritical behavior, this is not so as d* is greater than 3. To determine this 
dimension, we employ a simple scaling argument applicable whenever fluctuations 
dominate the critical behavior. We rescale the order parameter 77 by its bulk, 
uniform, value rjo and the position vector r by the correlation length £ in Eq.(6.26). 
Again, g and b are set to zero so that the Lifshitz tricritical point is approached. 
Requiring that the coefficient of the Laplacian squared term be independent of 
the reduced temperature t, we obtain r]o = £( 4_d )/ 2 ~ ^( d - 4 W 2 5 where d is the 
dimension of space. However ijq ~ t 13 by the definition of the exponent f3. Requiring 
the equality of these two exponents when the mean field values of (3 = 1/4 and 
v = 1/4 are inserted determines d* = 6. Below this dimension, fluctuation effects 
are important, and the exponents differ from their classical values. These exponents 
are as yet unknown theoretically as well as experimentally, so that a measurement 
of any of them would be of considerable interest. One reason why they have not 
been calculated is that the usual epsilon expansion for this system fails at dimension 
d=4. This can be seen by considering a term of the form rf in the integrand of 
Eq (6.26). Again rescaling r\ by 770 and lengths by £ one finds the coefficient of 
such a term to be £ d £p( 4-d )/ 2 . At d=4, every term of the form jf has the same 
coefficient £ d . This means that one should use the whole series in 77 to perform 
reliable calculations, which would result in an infinite number of renormalization 
group flow equations for the renormalized constants. 

Although the predictions of mean field theory fail at the Lifshitz tricritical point 
and within a region close to it, the predictions are valid outside this region, so that it 
is important to know its extent. This can be estimated with the help of the Ginzburg 
criterion. We follow de Gennes 3 ) and Joanny 5 ) and assume that fluctuations are 
unimportant if the fluctuations of the order parameter in a correlation volume is 
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much smaller than characteristic changes in the square of the order parameter itself, 
i.e. if 

< (S V ) 2 > /r,l oc l/r<% -> O)^ 2 « 1. (6.28) 
Therefore, mean field theory is valid for temperatures £* such that 

1 t du—r-2l3 

rf(^o)«~^/^ 

where ft, v and 7 are the mean field exponents. The N dependence follow from 
the fact that in polymers £ ~ iV 1 / 2 and T(0) ~ N, where N is the number of 
monomers in the polymer. Here we consider a system where all polymers have the 
same number of monomers. For d = 3, fluctuations can be ignored near the Lifshitz 
tricritical point provided that 

U > 1/iV 2 / 3 . (6.30) 

Similarly, 3> 1/N for an ordinary critical point (see section III) and t* 3> 1/iV 2 / 5 
for a Lifshitz critical point. In the limit N to infinity, mean field theory is valid in the 
entire critical region near the Lifshitz tricritical point. The influence of the Lifshitz 
on the surface tension of the coexisting A-rich and B-rich phases (see Fig. 23) is 
discussed in Appendix H. The elastic properties of the interface between the A-rich 
and B-rich phases can be find in Ref.[36]. 

VII. Rigid-flexible diblock copolymer system 
So far we have considered only the flexible systems. From Section II we have 
learned that also the rigid polymers can be described using this method. In this 
section we discuss the rigid-flexible diblock copolymer system. The order parameter 
densities have been already discussed in section II. D. They shall be used in the 
following subsections. 
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(6.29) 



A. The scattering intensity 

The scattering matrix has the following form in the case of rigid flexible diblock 
copolymers 50 -' : 

^(q)=(rf%) + rf> (q))~\ (7.1) 

Here 

(rf(q))- 1 = /(°)(q) = 



1 

V 



< 
< 



>A <PA >0 
>B 0A >0 



< 
< 



W >o 



< 
< 



A V a/3 >0 

(n)A(B) 

B V a/ 3 >0 



< 
< 



tin) MA) 

>A V Q/3 >0 

5 B V Q /3 >0 
)(B)A(A) 



<Q%W>o <QW>o <Q$'Q)2>o <Q\7Q" a 7>o 

< V a/3 <Pa >0 < V Q/3 <P B >0 < V Q/3 <^ 7 5 >0 < Q lS V Q/3 >0 



(7.2) 



is the ideal part of the inverse vertex function and is solely determined by the con- 
formations of a single copolymer molecule. Matrix (7.2) is a symmetric twelve by 
twelve matrix. By the direct calculations of the ideal averages, one may convince 
oneself that the averages involving the nematic tensor for the flexible chain, Q^(r), 
are N (N =N B or N^) times smaller than the other averages. This is understand- 
able since flexible chains consist of freely joined bonds, thus ordering one bond does 
not affect the ordering of the other bonds. In contrast, in the rigid rod, ordering one 
bond automatically orders the whole rigid B-part of the diblock copolymer. This 
allow us to neglect Q^fJ in our calculations and in the matrix(7.2). Thus 1^ is 
reduced to a seven by seven matrix. Futhermore, by chosing the coordinate frame 
in which the z axis is along the q vector, 1^ simplifies to the following block form: 



A 





A 



(7.3) 



Here 
A = 

1 

V 



'a <Pa >o 
'b <Pa >o 

] xx<l>A ^ ^0 

< Wyy 9a >o 



< 
< 



< 
< 
< Q 



A <Pb >0 

( n ) 2( n ) ^ 

B <I>B >0 
(B) ?(n) 
xx y b 



< <P A Vra >0 



< 



l(n)/)(B) 
'B ^ xx ^0 



< 
< 



>o <Qx B x ) Qx B x ) > <Qx B) Q[ B) 



l(n)A(B) 

) A ^ivv 

l(n)A(B) 

j b ^iyy 



>o 



'yy 



< Qyy 4>b ^0 ^ Qyy Qxx ^ Qyy Qyy ^0 



(7.4) 
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is a four by four matrix, whereas A is a diagonal three by three matrix, 





1 

V 



\ ^£xy ^xy -^0 











< Qxz ^ Qxz ^ >0 







(7.5) 



Because the off diagonal elements of Q v a/3 J appears only in A, they can be integrated 
out immediately in the integrals (see Appendix A,B). 

After the inversion of 1^ , we will make use of the incompressibility condition 
which, in Fourier space, reads 



^ n) (q)+^(q) = 0, 



(7.6) 



for q ^ 0. After these algebraic manipulations, we end up with the following form 
of the ideal and interaction vertex function: 



r (0) _ 

1 2 — 



A" 1 -2A" 1 +A 
A" 1 - A" 1 
A" 1 - A" 1 

iv 24 iv 14 



-1 A -l 
22 iv 23 



Ass' 
A34 



A2/ 



A3"/ 
A44 1 



(7.7) 



and 



■p(int) 



Po 



-2x 




4v B b -Ivbb 
-%vbb -3VBB 



(7.8) 



Here A~- is the ij element of the inverse of A. This matrix is also symmetric. 
To summarize this section, we have been able to reduce the set of twelve order 
parameters to the relevant set of only three order parameters, namely the two 

( B\ ( B} 

previously mentioned components of the nematic tensor i.e. Q xx (q), Qyy (q) and 
*(q). 

B. Stability limits of the isotropic phase 

At high temperatures the scattering matrix given by Eqs(7.1-8) is positive 
definite thus the system is stable with respect to compositional, and nematic per- 
turbations. However, as we lower the temperature the scattering matrix ceases to 
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be positive definite and the isotropic phase becomes unstable. We find the stability 
conditions by equating the determinant of the scattering matrix or its principal 
minors to zero. The results of this stability analysis are presented in Fig. 25 as a 
plot of xN versus the fraction of A monomers in the copolymer / = Na/N. If all 
phase transitions were continuous, these figures would be phase diagrams. First we 
discuss the case in which there is no explicit interaction tending to orient the rigid 
rods, vbb =0. This is shown in Fig. 25a. One can easily see that the nematic 
order parameters can be integrated out in this case, and the stability condition is 

An + A 22 + 2A 12 - 2 (A n A 22 - A 2 2 ) Po x = 0,. (7.9) 

The phase diagram encompasses in this case the disordered isotropic phase and 
some spatially ordered phases; for convenience we will call all of them lamellar 
phases. The stability analysis within our random phase approximation does not 
give us the symmetry of these phases. The stability condition against lamellar 
perturbations for the rigid- flexible system of the diblock copolymers (Eq(7.9)) has 
the same formal structure as that obtained in the Random Phase Approximation by 
Leibler 22 ) for the flexible-flexible diblock copolymer system. The difference between 
this condition and the one obtained by Leibler is in the ideal averages which, in our 
case, reflect the rigid- flexible structure of the polymer (see Appendix I). Let us 
discuss the differences in the stability limits for these two systems. First of all we 
note that in the rigid-flexible copolymer system the stability temperature is higher 
(or the interaction Flory-Huggins x parameter lower) in comparison to the systems 
of perfectly flexible system compare Fig. 25a and Fig. 12. For example below xN = 9, 
our system is stable for all f; while the flexible-flexible diblock copolymer system 
is only stable below = 10.5. This reflects the reduction in entropy due to the 
rigidity of part of the polymer. In addition, the mimimum of the curve shown in 
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Fig. 25a occurs at f=0.45 whereas in the flexible copolymers it occurs at f=0.5. The 
results of this calculations find support in the computer simulation 94 ) . 

For a system with nonzero anisotropic interaction parameter, i>bb, a nematic 
phase appears in addition to the lamellar and isotropic phases (Figs. 25 be)). The 
dashed line shown in Figs.25bc denotes the stability limits against nematic pertur- 
bations. It satisfies the following equation 

5 - = (7 10) 

2AT(l-/)2 2 1 } 

Furthermore we observe that the stability limit against lamellar perturbations goes 
towards higher temperatures as the anisotropic interaction parameter increases, 
although the change is rather small, a few percent at most. The black dot in 
the figures denotes the tricritical point at which the isotropic phase is unstable 
against both nematic and lamellar perturbations at the same temperature. As the 
transitions in the system are actually first order, this tricritical point is preempted 
by a triple point at which all three phases coexist. The calculated tricritical point 
is an estimate of the location of the actual triple point. 

C. The correlation functions 
We have calculated the following correlation functions for q=(0,0,g 2 ): 

GW(<fc) = y < *(qz)*(-qz) >5 (7.H) 

G* Q (q z ) = i < *(q z )Qi B A-q z ) >; (7-12) 

G QQ (q z ) = ^< $<?(&)$<!>(-&) >; (7.13) 
and the following Fourier transforms of these correlation functions: 
GW(*) = [ dx [ dy < (r)tf (0) > 

dq z G^(q z )exp(iq z z); 



1 

2tx 



(7.14) 
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G^ Q (z) = I dx I dy< tt(r)Q<f>(0) > 

■' J . (7.15) 

dq z Gm Q {q z ) exp(iq z z); 



1 

G QQ (z) = [dx [dy< Q<?(r)Q<f>(0) > 

1 /• ( 7 - 16 ) 

Here = — Q X x — Qyy^ measures the degree of orientational ordering of the rigid 
part of the copolymer along the z-axis. We note that as a direct consequence of the 
isotropy of the system the scattering intensities (Eqs(7. 11-13) can be rewritten in 
the following general form: 

Gw(q z ) = ^< tf(q*)tf(-q*) >; (7.17) 

G* Q (q z ) = ^ < *(q*)^Oq*)<$ >; (7-18) 

G QQ (q z ) = ^< q^fV^C^V)^ > ■ (7-19) 

Here q* is an arbitrary vector of length |q*| = |g^|; q* is a unit vector along q*. The 
summation over repeated indices is implied. 

The calculations have been performed in the limit q z — > 0, Na,N b — > oo, 
such that q z N B =const and / = Na/(Na + N B ) =const, and for the following 
sets of parameters: (1) f=0.20, %iV =19.51, v B b/x =0-4; (2) f=0.20, %iV =19.51, 
vbb/x = 0; (3) f=0.62, X N =11.35, v BB / X =3.0; (4) f=0.62, X N =11.75, v BB / X = 
0. In all cases the parameters has been chosen such that the system is close to the 
stability limits which are: for case (1) xiV=19.53; for case (2) x-^=19.56; for case 
(3) xN =11.40; and for case (4) X N =11.79 (see Fig.25). 

In Fig. 26 GW((/ 2 ) and G^^{z) are shown for f=0.20 and f=0.62. In the case 
of a long rigid part, f=0.2 (Fig. (26a)), there are two peaks at q\ and qi, whose 
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positions are determined by the length of that part i.e. q\ 2tt/((1 — f)Nl) and 
q 2 ~ 47r/((l — f)Nl). In the case of a short rigid part, f=0.62, there is only one peak 
at qi (Fig(26c)). The presence of the second harmonic in Fig(26a) and its absence in 
Fig(26c) indicates that the domain boundaries are sharper in the system containing 
longer rigid pieces. This is supported by the plots of the real space correlation 
functions. For f=0.62 (Fig.(26d)) this function is almost sinusoidal, whereas it has 
additional structure for f=0.20 (Fig. (26b)) In both cases the length of the rigid part 
of a copolymer sets the length scale. Finally by comparing the cases of vbb/x 
(solid line in Fig. (26)) and vbb/x = (dashed line in Fig. (26)) we note that the 
general features of these functions are not affected very much by the Maier-Saupe 
parameter vbb- 

The nematic structure function, (jqq(<? 2 ), and its Fourier transform (the ne- 
matic order parameter — nematic order parameter correlation function), Gqq(z), 
are shown in Fig(27) for f=0.62. The peak at q z = reflects the tendency toward 
nematic order, while that at q z ~ 2tt/ fNl reflects the tendency to lamellar order. 
They appear in Gqq(z) as oscillations due to the nascent lamellar order superim- 
posed on the decaying nematic order. Not surprisingly, the Maier-Saupe interaction 
strongly affects the nematic signature of this correlation function as comparison 
with the case vbb/x = clearly shows. For f=0.2 the tendency to nematic order is 
greatly increased (because the system is closer to the isotropic-nematic instability) 
so that GQQ(q z ) is completely dominated by the peak at zero wavevector, and the 
oscillations due to nascent lamellar order are ignorable on the same scale. For this 
reason, we have not shown this function or its Fourier transform. 

In Fig. (28) the mixed correlation functions G^Q(q z ) and G^q(z) are shown for 
f=0.20 and f=0.62. Just as for the density-density correlation function (Fig. (26)), 
the structure of the density-nematic order parameter correlations is richer for f=0.20 
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than for f=0.62. Similarly, these structures do not depend very strongly on the 
Maier-Saupe parameter. The most interesting observation to be made comes from 
a comparison of the functions Gw(z) and G^,q(z) for f=0.62 (Figs.(26d,28d)). We 
see that the nematic order parameter in a given direction and the linear concentra- 
tion of rigid rods in that direction are almost exactly anticorrelated; i.e. large and 
positive Qlz (z) occurs in planes normal to z which are rich in the flexible A-part 
of the copolymer. At first sight, this result may appear to be counter intuitive. 
One knows, for example, that in a mixture of rigid and flexible homopolymers, a 
higher density of the rigid polymers is found in the nematic phase. However we are 
concerned with the isotropic phase, and the effect can be understood as arising from 
that isotropy. As an example, let us assume that we are at a point in space in which 
the density of B monomers is greater than that of the A. For simplicity, we will as- 
sume that there are B monomers only out to a distance R, and A monomers from 
R to bR, with b = (1 — /) -1 / 3 , so that the relative concentration of A monomers is 
f. With such a density distribution, 



~ / ds27r S P 2 (-==)).(7.20) 
Carrying out the integrations and setting R = NbI/2 = N(l — /)//2, we obtain 



9r l , , 2x ,,3, 1 



PolGv Q (z) = F(z/N B l), 

F & = 18x2[ l - ln( (T^7)^ )] + ?(TF " 2) - (7 - 21) 

As long as the fraction of flexible monomers is not too large, / < 1 — 2~ 3 / 2 ps 0.65, 
this function is negative at zero argument, thus reproducing the anticorrelation of 
density and nematic fluctuations. Further, it has a maximum when the argument 
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is approximately 0.36(1 — f) 1 ^ 3 . In the above picture, we expect the spacing of 
spheres, and hence the period of the correlation function, to be N B l/(l — f) 1 ^ 3 . As 
the model correlation function is minimum at z = 0, it should also be minimum 
at integer multiples of this period with maxima a distance 0.36(1 - f) 1 ' 3 on either 
side. There are additional minima at half odd integer multiples of the period. For 
f=0.2 as in Fig. 28b, this model predicts minima at z/lN B = 0, 0.54, 1.08, 1.62, 
2.16, etc., and maxima at z/lN B = 0.33, 0.75, 1.41, 1.81, etc. These numbers are 
in satisfactory agreement with the figure. The double peak structure which this 
model, with its sharp density profile, produces requires at least the presence of 
a second harmonic to approximate it. As we have seen, the system with shorter 
rigid sections, f=0.62, does not exhibit a second harmonic in its structure function. 
Thus the above model does not describe the function of Fig. 28d very well at all 
except to reproduce its periodicity and the important negative value at the origin. 
The above model should not be emphasized unduly, but it does serve to emphasize 
that the anticorrelation of density and nematic fluctuations that we have found is 
reasonable. 

The behavior of all correlation functions close to tricritical points is similar for 
all f>0.25, and in particular to the case f=0.62 presented here in detail. For fc0.25 
the correlations functions begin to develop the structure shown here for f=0.2. 
Finally we note that the locations of the minima and maxima of the correlation 
functions are relatively independent on the Flory-Huggins parameter. Changing it 
changes the heights of the peaks, but hardly affects their positions. Thus the main 
features of the local domain structure are determined by the geometry (rigidity and 
length) of the copolymer molecule and excluded volume effects. In our calculation, 
the latter is modelled by assuming incompressibility of the system. 
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VIII. Rigid-flexible blends 



Here we shall show that within the approach presented in this paper we can 
also obtain the general features of the nematic phase in the system containing rigid 
and flexible main-chain nematogenic polymers. The other systems like the mixtures 
of two different rigid polymers, and a system of n-block copolymer consisting of n 
rigid and flexible parts are described elsewhere 60 ) . The Landau- Ginzburg expansion 
of the free energy is derived, and we show that the general feature of the isotropic 
nematic phase transition in the system of short flexible polymers and long rigid poly- 
mers are qualitatively the same as obtained in the completely different treatment 
of the density functional theory 95 ' 96 ) applied to infinitely long rods in the solvent. 
The latter model is commonly known as the Onsager model 97 ) The detailed study 
of the (extremely rich) phase diagrams produced in the Landau-Ginzburg model 
and its variants is contained in Ref [57, 98-100]. 

In order to describe the system we shall use two order parameters: the usual 

( B) 

^ parameter and a nematic order parameter Q B = Qzz • The nematic order 
parameter for the flexible part, Sa = Qzz , can be neglected for the reasons already 
stated in the previous section and in Ref [60] . 

The Landau-Ginzburg free energy for the rigid flexible mixture (with the frac- 
tion of the monomers in the flexible polymers being 0, can be easily calculated using 
the approach presented in Appendix A and B. We find the following form of the 
free energy per monomer: 



(8.1) 



0/ is given by Eq(3.5); 



^b(Qb) = 



5 



vbb 
2 



Qe 



25 



- 2 Q% + 



425 



Qb- 



2N B (l-</>) 



21N B (l-4>) 



196A^ B (1 -0) 3 



(8.2) 
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is the part of the free energy coming form the nematic ordering of the rigid polymers; 

(8.3) 

is the term, representing the coupling between the order parameters in the mixture. 

Given the free energy Q(Q B ,^) : one easily finds the conditions for thermal 
equilibrium of two or more coexisting phases. The condition 

= (8.4) 

oQb 

must be satisfied in the absence of external field in all coexisting phases. Addition- 
ally the function 

G(Q B , *) ee Q(Qb, *) - ^§§^ *- (8-5) 

must be equal in all phases. 

For the two phase coexistance the equations are solved for the three unknowns, 
Q B \ Qb , an d ^ 5 i- e - the values of the order parameters in the two phases. In the 
isotropic phase Q B = 0, but as we shall we we have the cases when two nematic 
phases of different Qb can coexist. Please note that at coexistance we have in one 
phase * = 0. In the other phase \1/ is equal to the difference in the composition of 
the coexisting phases. 

We consider the mixture of flexible A and rigid B polymers in which the latter 
are much longer than the former, Nb Na- At high temperatures one finds that 
the composition difference between the B-rich nematic phase and A-rich isotropic 
phase is extremely small, namely 

4>n - h ~ ^ « 1, (8.6) 
1Mb 

where 0/ (4>n) is the concentration of A monomers in the isotropic (nematic) phase 
at coexistence. Neglecting terms of the order of Na/Nb or smaller, the two-phase 
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coexistence region shrinks to a line given by 

l«O.21vj^(l-0)JV B . (8-7) 

where v B N B is the interaction parameter at the isotropic nematic phase transition. 
Please note that the parameter is rescaled by ksT. The nematic order parameter 
jumps from zero in the isotropic phase to Qb ~ 0.275(1 — 0) in the nematic phase. 
Here <p is the average concentration of A monomers in the mixture at N-I coexistence. 
Recalling that Qb is the nematic order parameter of the B monomers divided by 
the number of all monomers, we see from this result that the jump in nematic order 
parameter per B monomer is constant, independent of the concentration of the A 
component. These results are similar in spirit to the results obtained by Warner and 
Flory 101 ) and Onsager 97 ) for the system of rods diluted by a small particle solvent. 
If we assume that the monomers in the rigid polymer interact only via excluded 
volume interactions we can write vbb ~ PqVq and po(l — 0) = PbNb, where ps is 
the density of rigid polymers and vq is proportional to the excluded volume of the 
two monomers. Since vq ~ / 2 -Ds, where / is the length of the monomer, Lb = INb 
is the length of the rigid polymer and Db is the thickness of the polymer we finally 
find the famous Onsager result that in the limit of the rigid rod much larger than 
the solvent (either monomeric or polymeric) and Lb ^> -Ds, the rod density at the 
isotropic-nematic transition is inversly proportional to L 2 B i.e.: 

p B L%D B ~ 1 (8.8) 

Since the Onsager scaling is exact, Eqs(8.7-8) prove that the method applied here 
to rigid polymers is correct at the fundamental level. 

We also find the phase separation inside the nematic phase. This phase sep- 
aration is goverened by the Flory Huggins interaction parameter. The location of 
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the critical point is given in the limit of N B 3> Na by (see Eq(3.6,3.7): 



mf 



c 



2N A ' 



1 



4>? f = i- 




(8.9) 



As is evident from Eq.(8.9) the nematic ordering does not change the critical point 
temperature, which is the same as in the ordinary polymer blend discussed in section 



In this paper we have applied the Landau- Ginzburg model to several systems: 
the polymer blends, diblock and multiblock copolymers, rigid and flexible polymers 
and their mixtures. The Landau- Ginzburg free energies and scattering intensities 
were calculated from the Edwards mesoscopic hamiltonian using the field theory 
methods. We hope that this review will be helpful for the future studies of polymer 
melts and blends. 
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Appendix A: The cumulant expansion 

The method of cumulant expansion is applied here to the binary mixture of 
A,B homopolymers, but can be easily generalized to any other polymer system. 
The partition function, Z\(^^\ as a functional of prescribed Fourier transform of 
the concentration distribution, (f)y \<i), (7 = A, B)is equal to 

Z[<p^]=exp(-H[ ( p^]/k B T)( J] 5(0^)(q)-^)(q)} o (A.l) 

1=A,B 

where the average over the conformations of the chains is 

< • • • > = J\ Jj / Dv- Dv 1 ■ ■ ■ W A [r a ]W B [r p }. (A.2) 

a=l p=l J 

and the interaction hamiltonian, H, (see section II) is 



H[^\ =k B T P0 J ^^aa\^ a) {^ 

+\^bb\^ b) ^)\ 2 + w AB ^ A) (q)4 ns) (-q))) • 



(A.3) 



Now, we can introduce the integral representation of the delta function in 
Eq.(A.2) 

n^(^ ) (q)-^ ) (q)) = 

J] / DJ,e W (i | (^>(q) - #^>(q))J 7 (-q)) 

and expand the exponential in cj)^\ Eq.(A.l) now reads 

z[^"-')]=exp(-H / [^"'' ) ]/*Br) II (/w 7 ex P ^|^3^>(q)^(-q) 



X 

n=0 



(A.5) 
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where Sl^ denotes the n-th moment 



'n 



^ (7) (qi"'-qn) =(^ ) (qi)---^ ) (qn)) Q . (A.6) 



It is obvious that ideal averages over A chain and over B chain can be performed 
independently, thus there are no mixed moments. The mixed moments are found 
only if the chains are connected inside a molecule. This is the case of the diblock 
copolymer system. 

The logarithm of the moments expansion is the cumulant expansion 7 > 14 ). Thus 
the second term in Eq.(A.5) can be rewritten as 

J] J DJ 7 ex P (F^[jj + i J (03^ } (q)^(-q)), (A.7) 

where the cumulant expansion, F 7 n7 ^[J 7 ], which is a functional of the coupling 
fields, J 7 , is given by 

^ \M-^, n[ J (2yr)3 J {2n) 3 b n (qi qnj 
n=l 

5(qH — + qn)^ 7 (-qi) • • • J 7 (-q n )- 

It is easy to represent cumulants in terms of the moments 14 ' 22 e.g. 

*(qi + q 2 )4 A) (qi,q 2 ) = ^2 (A) (qi,q 2 ) - 5 1 * ( ^ ) (q 1 )^ ( ^ ) (q 2 ). 
It is also easy to compute the moments and cumulants explicitely, e.g. 

i<^(q)^(-q) > = 
tiaNa ( 1 \ f (I + a a 



(A.8) 



PoV 

where a a = sinqlA/qlA and q 7^ 0. Expanding this expression in qlA we find 

1 < <^(q)0^(-q) > = ^^g(x A ) + N A 0(1/N A ), (A.9) 
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where = uaNa/poV is the average concentration of A monomers in the system, 
xa = Na^Ia) 2 /6 and <7(x) = 2(x + exp(— x) — l)/x 2 is the Debye function. Since 
we are interested in 1/ yjN 1 corrections we can neglect all higher order corrections 
(e.g.l/iVA corrections in Eq(A.9)). To this point, no approximations have been made 
other than that which enabled us to express the conditional partition function in 
the form of Eq.(A.l) 

Appendix B: The saddle point approximation and RPA 
Now the integral in expression (A. 7) is calculated by the method of the steepest 
decent, i.e. we approximate (A. 7) as follows. 



J Dj 7 ex P (V(^)[j 7 ] + z J (03^(q)^(-q)) « 

exp (F^[j;] + i J ^^>(q)J*(-q)) 
Here, at the saddle point, the coupling field, J*(q), satisfies the equation 



(B.l) 



(n 7 ) f 

t^(q). (B.2) 



Thus, to within an uninteresting constant, the expression (A. 7) is approximated by 
J] exp(-F^>[* 7 ])= J] exp(^[J 7 *]+* /^^(^(-q)) . 

~1=A,B J=A,B ^ J \ ) J 

(B.3) 

Note that Eq.(B.3) together with Eq.(B.2) constitute the Lengendre transform. We 
have F^ } [tf 7 ] : 

X(5(qi + • • • + q n )* 7 (-qi) • • • * 7 (-q„), 
with the coefficients, which are simply expressed in terms of the cumulants 7 ' 14 ' 22 
e.g. 

r 2 A) ( qi ,q 2 ) = l/^ A) ( qi ,q 2 ). (B.5) 
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(B.4) 



A j (qi,q 2 ,q 3 ) = --^ ^j- (B.6) 

s 2 (qi,q2)^ 2 (q2,qs)^ 2 (qi^qs) 

The same formulas are obtained for B-chains. The formulas for functions for 
n > 3 are rather lengthy. Here we would like to give the abbreviated formulas for 
n = 4, 5, needed in our later calculations. They are as follows: 

r< A > = - (s™ - 3 J s< A >r< A >s< A >) (r^) 4 , (B.7) 
4 A) = - (s {A) + 15 J 4 A) rf } / s^4 A) s { 3 A) -io J s[ A) r^s {A ^ (ri A) ) 5 , 



(B.8) 

where the proper symmetrization and the integration over dummy variables are 
assumed. 

We note that Eq.(A.l) together with Eq.(B.4) constitute the basis of the ran- 
dom phase approximation (RPA) 2 ' 22 ) and the RPA vertex functions are easily 
obtained by summing the interaction term form the hamiltonian and the term in- 
volving and . For example we have 

rf } = + rf> - 2 XP0 (B.9) 

and 

r (0) =r (A)_ r (B)_ (B1Q) 

In the limit of q — > we find 

r^ ) (o,o,o) = -ir^ ) (o,o) (b.ii) 

(analogously for B with changed to 1 — 0). The right hand side of Eq(B.ll) is 
given by the inverse of Eq(A.9) (see Eq(3.9)). 
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Appendix C: One-loop equations for the vertex functions 

Fig. 7 and Fig. 8 show the one-loop diagrams contributing to T 2 and r 3 , re- 
spectively. Note that vf ] and vf ] contribute to T 2 (Fig. 7), while r^ 0) and r^ 0) 
contribute to r 3 (Fig.8). The solid lines represent l/T^ while the higher order ver- 
tex functions are represented by points in the diagrams. The n-body vertex function 
has n lines emanating from the point representing this function. The dashed lines 
are not parts of the diagrams but are included only for clarity (according to Ref.[7]). 
The final approximation (Hartree) in the framework of the loop expansion is to make 
the equations self consistent. Here it corresponds to the change of r 2 and Ig 
appearing in the integrals, represented by the diagrams in Figs. (7,8), to T 2 and T 3 . 
The final self-consistent one- loop equations for the two-body and three-body vertex 

functions are as follows: 

t~i / x -p(O), x . 1 f dk rf^q, -q,k, -k) 
T 2 (q, -q) = F 2 >(q, -q) + - J ^ ^ _ k) 

1 f dk (r 3 (q,k,-q-k)) 2 



(C.l) 



2 J (27r) 3 r 2 (k,-k)r 2 (q + k,-q-k)' 

r 3 (q,p, -q-p) = r^ 0) (q,p, -q- p) 

dk r 3 (q, k, -q - k)r 3 (p, k, -p - k)r 3 (-q - p, p + k, q - k) 



+ 



(2tt)3 r 2 (k, -k)r 2 (q + k, -q - k)r 2 ( P - k, -p + k) (C.2) 

i(0), 



+ 1 f dk r^(k,-k,q,p,-q-p) 



2 J (2tt)3 T 2 (k,-k) 
Before proceeding, we first want to determine the order (in terms of A^ 7 , 7 = A, B) 

of successive terms in Eqs.(C.l,2). For simplicity, we consider for a moment the 

case of Na = Nb = N and I a = Ib = I- For equal length of A and B polymers 

the critical concentration is C = 1/2 and at the critical point only the first and the 

second terms in Eq.(C.3) are nonzero. By comparing these two terms we find the 

following relation at the critical point: 

1 r dk rj 0) (k,-k,o,o) _£_ . . 

\Tf (0,0)1 J W T 2 (k,-k) ~ A + A3' [ ■ > 
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where 2n/A is the upper cutoff in the integral. As we can see the choice of the 
cutoff is crucial for determining the relative magnitude of the mean-field value of Y2 
(that is r^ ) and the first (one-loop) correction to it. There are no ready recipes 
for the choice of the cutoff. Here we take it proportional to the radius of gyration 
i.e. to VW. The discussion of the cutoff is given in Appendix D. 

In the solution of Eq.(C.l,2) we use the following approximation for Ij (sim- 
ilarly for T 2 ') (see Appendix A and B): 

r^( q ,- q) = i±Jp. (c.4) 

(pN A 

The coefficient (Eq.C4) in front of q 2 is obtained from the expansion of the expres- 
sion (A. 9) for large qlA- Finally we observe that 

ifV -q,0) = -r^(q, -q)/0, 

rf } (q, -q,0,0)=2r^(q, -q)/0 2 , 

4 A) (q, -q, 0, 0, 0) = -6rf } (q, -q)/0 3 . 

We have analogous equations for B chains with changed to 1 — 0. 

Appendix D: The upper wavevector cutoff in polymer blends 

As we have seen in section III the fluctuation corrections to the critical tem- 
perature strongly depend on the upper wavevector cutoff. Although it is ubiquitous 
in statistical mechanics there are no ready recipes for the choice of the cutoff. In 
low molecular mass liquids we usually have one natural length scale, which corre- 
sponds to the size of a molecule and the cutoff is usually made proportional to this 
length scale. In the polymer mixtures the problem is more complicated since we 
have three different length scales: the total length of a polymer molecule, Nl, the 
size of the region occupied by a polymer molecule (proportional to the radius of 
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(C.5) 

(C.6) 
(C.7) 



gyration), ~ y/~Nl, and finally, the microscopic length scale /, which is determined 
by the size of a single monomer. We believe that A should be proportional to the 
radius of gyration, that is to y/~N. At high temperatures, where the interactions 
between the monomers are irrelevant, the monomer- monomer correlation function 
decays exponentially with the characteristic correlation length which is proportional 
to the radius of gyration. Moreover, in all our calculations the specific structure 
of monomers has not been taken into account. The microscopic length scale is ir- 
relevant here. Finally, if we chose the microscopic length scale as the cutoff the 
fluctuation corrections would survive in the limit of N — > oo and thus RPA would 
not be the correct description in this limit. We believe this is not the case. We 
note that since the size of the polymer molecule in the blend is roughly propor- 
tional to ; this choice is also in accordance with the prescription known from 
low molecular mass systems where the cutoff is made proportional to the size of a 
molecule. In the case of assymetric mixture we postulate that the cutoff should be 
equal to 2tt / (A(Ra, Rb)), where A is the symmetric function of the two radiuses of 
gyrations for A and B chains. Finally, we assume the scaling form for A : 

A(R A ,R B ) = R A f(x), (D.l) 

where x = Rb/Ra- From the symmetry properties of A, we have the following 
equation for f(x): 

§/<»)-/(;)• (D.2) 

The general solution of this functional equation is 

f(x) = ^C 1 1 -±^ + C 2 x. (D.3) 

Here C\ and C<i are two constants which depend on the microscopic details of the 
system. We cannot determine them from the present theory, thus, in the interpre- 
tation of experimental results, they should be the fitting parameters. 
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In the simplest approximation we can set Ra ~ VWaIa- In general in the 
polymer system this radius changes with temperature and concentration and so 
it would be desirable here to determine it self consistently from the theory. The 
equation for the radius of gyration, R, offers a simple way for the self consistent 
determination of the cutoff, if we set it equal to 2nC/R where C is some numeric 
constant. Thus apart from this constant the cutoff would be determined from the 
theory. In section III D the equations for the radius of gyration are shown to involve 
the collective structure factor (proportional to the scattering intensity Eq(3.9)). In 
this way the equation for the collective structure factor and the equation for the 
radius of gyration are mutually coupled. 

Appendix E 

Here we present the results used in section IIIC for the calculation of the radius 
of gyration. Here we shall follow the equations given in Ref [20] . For the fourth order 
correlation function we have: 

< jff (q) iff (-q) iff (k) iff (-k) > = s 9l (« , o, k) + 2( 9l ( q , q + k, q ) 

+ gi (q, q - k, q) + 9l (k, q + k, k) + gi (k, q - k, k)) (E.l) 
+ %i (q, q + k, k) + g x {q, q - k, k)). 



Here 



N A p(xi) 
9i (gi, 92, qs) = + 



XxX 2 Xs x 1 (x 1 - x 2 ){xi - x 3 ) 
p(x 2 ) p{x 3 ) 



(E.2) 



x 2 (x 2 - x 3 )(x 2 - xi) x 3 (x 3 - x 2 )(x 3 - Xi) ' 
where Xi = qfl 2 /6 and p(x) = (exp (— xNa) — l)/x. For the two point correlation 
function we find 

<02 ) (q)02 ) (-q) >o=2(Na+p(x))/ Xi (E.3) 

and x = q 2 l 2 /6. 
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In order to obtain the radius of gyration one has to differentiate 
fc(k,q) =< ^\^\- q )^\k)^\-k) > 

(E.4) 

- < ^(^(-q) >o< ^(k)^(-k) > 
twice with respect to q and finally take the limit of q — > 0. We find 

/i"(k,0) =/(k) = 4Nl cos 2 6 (-120 + 120 exp (y) - 96y - 24 exp (y)y - 24y 2 
-12 exp (y)^ 2 + 4 exp (y)?, 3 + y 4 )/(3 exp (y)?, 5 ), 

(E.5) 

where y = NAk 2 l 2 /6 and 6* is the angle between the q and k vector. 

Appendix F: One-loop equation for 1?2 in the lamellar phase 

In the lamellar phase of the diblock copolymer system the field $ is approxi- 
mated by the single harmonic: 

* = Acos(qiz) (F.l) 

In general qi, which describes the periodicity of the phase, is different from the q* , 
which marks the divergence of the structure factor. In practice they are very close. 
The vertex function in this case, T 2 (q, A) is the function of the amplitude A and the 
wavevector q. In the disordered phase A=0, and T 2 depends only on the modulus 
of q. The equation for the vertex function for the symmetric case (/ = 1/2) is as 
follows 20 ^: 

r 2 (c l ,A) = rf\\c l \,A = o)+ 1 - J ^ I ri°>(q,- q ,k,-k)rj 1 (k,A) 

+ A 2 rf\c l ,-c l ,q l ,-q l ) 

Please note that the function, T 2 is not isotropic in the lamellar phase, because a 
particular direction (z in this case (Eq(F.l))) has been selected. The dependence 
of on angles is rather weak. 
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Appendix G 

The coefficent in the Landau- Ginzburg free energy Eqs(6. 26-27) are as follows: 

a = 71 

, 27273 - 74 



4 73 

473 75 - 3737478 + 3737I76 - 7477 

127! 



(G.l) 



where 



7i 
72 

73 

74 

75 
76 
77 



1 X N 



(1-0) 5 
1 

3(1 -0) 3 ' 
1 

0(1-0)' 
1 

Mi 1 ' 
1 

5(1 -</>) 5 ' 
1 

(1-20) 



78 = TT—nl- (G.2) 



2(1 - (f)) 2 (j) 2 
1 

(1- 

Note that with the exception of 71, all coefficients are determined solely by the 
entropy of mixing. 

The consolute line is given by the vanishing of the coefficient a in Eq. or, from 
Eqs. (G.l) and (G.2), 

* N = jr*y (&3) 

A tricritical point arises when the coefficients a and b vanish simultaneously, which 
occurs at the point 

XtriN = 6, 
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<t>tn = 2/3. (G.4) 

We have established, therefore, that this point, the point at which the disorder line 
and Lifshitz lines all meet, is a trier itical point. To establish that it is a Lifshitz 
tricritical point we have to check go. Here we have 

90 = -Ajt^T)-^^)—> (a5) 



1/1 9( X N) 2 <P ( X N) 3 <P 2 \ iV 2 Z 4 
= 36V(T^T 16 8 J"36— (G - 6) 



We immediately verify that at the tricritical point Nxtri = 6, 4>tri = 2/3 at which 
the coefficients a(x, 0) and &(%, <fi) both vanish, the coefficient go(x 0) a l so vanishes. 
Therefore the point is a Lifshitz tricritical point as stated (see also Appendix H). 
Appendix H: Influence of Lifshitz Tricritical Point 
At values of x larger than those given by the consolute line (in section VI), 
the system separates into A-rich and B-rich phases (Fig. 23). We shall calculate 
approximately the surface energy between these two phases to illustrate the effect 
of the proximity of the Lifshitz tricritical point. The two coexisting phases are 
characterized by uniform order parameters of amplitude ±770, where this amplitude 
is obtained by minimizing the free energy of Eq. (6.26) and requiring a uniform 
solution. One obtains 



—b + Vb 2 — 3ac 



ni = Fc • (H.i) 

To determine the interfacial profile rj(z) between these phases, we require a non- 
uniform solution of the Euler-Lagrange equation which expresses the minimization 
of the free energy Eq. (6.28). The solution must also satisfy the boundary conditions 



rj(±oo) = ±T] . 
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(H.2) 



The Euler-Lagrange equation is of fourth order and nonlinear, and we are unable 
to solve it. In order to obtain an analytic form for the surface free energy, we 
approximate the bulk part of the free energy by a double parabola with the same 
curvature at the minima; 

-( \a\r, (r) + br, (r) + cr? (r)j ~ j 2 \a\y(s)( V (r) + Vo ) 2 if 77 < 0, (H ' 3) 

where 

y(s )= [( l + ,2)1/2 _ s](l + s2) l/ 2) (H _ 4) 

and 

s = b/(3\a\c) 1/2 . (H.5) 

The parameter s is a scaling field, a dimensionless measure of the deviation, 6, 
of the system from tricritical behavior measured with respect to a 1 / 2 , where a is 
proportional to t = (T — T con )/T con , the deviation of the temperature from the 
consolute temperature. With this approximation in the free energy, the Euler- 
Lagrange equation is piecewise linear 

ho d^~ 9o d^ + 4 ^ y ^ ± ^ = a (H - 6) 

The form of the solution depends upon the ratio y(s)/p 2 where p is a second scaling 
field 

p = W(4|aM 1/2 , (H.7) 

a dimensionless measure of the amplitude of the square gradient term in the free 
energy. It is a measure of the deviation of the system from Lifshitz behavior. The 
ratio of y/p 2 is easy to interpret; the scaling field p for fixed go diverges as the 
consolute line is approached. Similarly, the scaling field s for fixed b also diverges as 
the consolute line is approached. However, the function y(s), only varies by a factor 
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of two as s varies from zero to infinity. Thus y/p 2 will always be smaller than unity 
sufficiently close to the consolute line; how close in terms of temperature depends 
on the amplitude go, that is, how close the system is to Lifshitz behavior. 
In general the form of the profile is 

n ( z ) = { -Vo + K 1 exp(z/£,+) + K 2 exp(z /£,-), if z < 0; , . 

Iy ' \ T] -K 1 exp(-z/Z+) - K 2 exp(-z/£-), if z > 0. 1 



with 



Kl = (e + /e-) 2 -i' (IL9) 
* 2 = - (wfc' (H - 10) 



and two decay length in the system 

1/2 

h J r v--^- p 



i 2 = 2 m p (i± (i-^). ( H. n) 



As the consolute line is approached, the correlation length £_ diverges as t -1 / 2 , while 
is finite unless the system approaches the Lifshitz tricritical (or critical) point. 
It is easy to see that for y/p 2 < 1, £± are real and the profiles are monotonically 
decaying functions at large distances from the interface. For y/p 2 > 1, that is, 
sufficiently close to the Lifshitz tricritical point, £± are complex and the profiles 
are no longer monotonically decaying, but are damped oscillatory functions. Using 
Eq(H.8-H.ll) this oscillatory profile can be written in the following form: 

, v f -7/o +exp(z/f)(77oOOs(27rz/A) + K 3 sin(27T2;/A)) if z < 0, ,„ * 
l{) \vo+exp(-z/0(-Vocos(2nz/\) + K 3 sm(27iz/\)) if z > 0, 1 J 



where 



* = 7s^ <H ' 13) 
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and 

Note that this profile has the similar form as the correlation function of the disor- 
dered phase on the copolymer-rich side of the disorder line (section VI B). We note 
here that the non-monotonic profiles are encountered in homopolymer-rich phases 
sufficiently close to the Lifshitz tricritical point. The correlation length £ diverges 
as the Lifshitz tricritical point is approached as t -1 / 4 , as expected from Sec. VI C. 

The interfacial tension is obtained by inserting these profiles into Eq. (6.26) 
and dividing by the area of the planar interface. The result can be written in the 
scaling form 



where 



y 



3/2 (2P + V 1 / 2 ) 



y (q n) _ J (1 + S 2 ) 1 / 2 [ P ~ (P 2 - y)V2]l/2 + [p + (p2 _ y )l/2]l/2 
(l + s 2)l/2 ( p + ?/ l/2)l/2 



if y < P 2 

if y > p 2 - 

(H.17) 



The surface tension is proportional to iV -1 / 2 as ho is proportional to iV 2 . This 
expression for the surface tension has the following properties; on approaching the 
Lifshitz tricritical point, t — > 0, p and s constant, 

a(M,p)~t 5 / 4 ; (H.18) 

on approaching an ordinary tricritical point, t — > 0, s constant, p ~ t -1 / 2 , 

a(t,s,p)~t; (H.19) 

on approaching a Lifshitz critical point, t — > 0, s ~ t -1 / 2 , p constant, 

a(t, S ,p)~t 7 / 4 ; (H.20) 
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on approaching an ordinary critical point, t — > 0, s/p constant, s ~ t 1 ^ 2 , 

<r(t,8,p) ~ t 3/2 . (H.21) 

All of these behaviors are the ones expected in mean field theory. The crossover 
between these behaviors is governed by the two scaling fields s and p of Eqs. (H.5) 
and (H.7). 

Sufficiently near a transition so that fluctuations become important, the surface 
tension will still have a scaling form 

-(^o)=^(^,^), (H.22) 

but the form of the scaling function E will differ from that given above, and the 
critical exponent \i and the crossover exponents 4>t and 4> g will differ from their 
mean-field values 5/4, 2, and 2, respectively. As noted earlier, there are no calcula- 
tions as to the values these exponents should take. 

The presence of the Lifshitz tricritical point also affects the results of scat- 
tering measurements as seen in the structure function. At long wavelengths, and 
sufficiently close to the consolute line, the form of a general structure factor is 

I(q) = (a + g q 2 + Kq 4 )-\ (H.23) 

where K, a function of the copolymer concentration, will depend on the particular 
structure function. On approaching any point on the consolute line except the 
Lifshitz tricritical point, it is convenient to rewrite this expression in the form 

r( \ -1/1 , 9o 2 . Ka 9o 4\-l 
7(g) = a (H q + — — ^q 

a g£ a 2 

= a- 1 (l + (q0 2 + K'p- 2 (q0 4 )-\ (H.24) 

where, to within numerical factors, £ is the ordinary correlation length £_ of Eq. 
(H.9). As the scaling field p diverges as the consolute line is approached, the term 
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in g 4 becomes unimportant and the structure function can be written in the scaling 
form 

I(t,q) = ts(q£) 
s(x) = (1 + x 2 )- 1 , p^oo (H.25) 

that is, the usual Lorentzian form. On the other hand, if one approaches the Lifshitz 
tricritical point so that the scaling field p remains finite, it is more convenient to 
write Eq. (H.23) as 

= a" 1 (l +A-"p« 9 ) 2 + (5,) 4 )- 1 , (H.26) 

where K" is an unimportant function of the copolymer concentration and £ is again 
the correlation function which behaves as t -1 / 4 as the Lifshitz tricritical point is 
approached. Again, this can be expressed in a scaling form 

I(t,q) = ts(qO 
s( x ) = (1 + K"px 2 + x 4 )-\ (H.27) 

with p fixed. Therefore, the presence of the Lifshitz tricritical behavior makes itself 
known by the fact that, in a fit of the structure function at long wavelengths, terms 
of order q 4 cannot be neglected with respect to terms of order q 2 . Again, when 
the effects of fluctuations are included, one expects the structure function to take 
a scaling form 

I(q, T, 6, g ) = PsiQtfr ^ 7^ (H.28) 

but the form of the scaling function will differ, and the critical exponent 7 and 
the crossover exponents <pt and <p g will differ from their mean-field values of 1, 
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2, and 2. Experiments on mixtures of homopolymers have observed the change 
in the exponent 7 from its mean-field value of 1 to its fluctuation driven value 
of 1.26 as the consolute line was approached (section III). Presumably, a similar 
change could be observed if the Lifshitz tricritical point were approached, because 
the region of fluctuation dominated behavior is somewhat larger. Again, such an 
observation would be most interesting as there exist no calculated value for this 
exponent characterizing the Lifshitz tricritical point. 

Appendix I: The ideal averages for the rigid-flexible diblocks 
Below we show the results of the calculations of the ideal averages. Here a = 
sin(ql)/ql, f is the fraction of A monomers in the copolymer and N=N^+N B is the 
number of all the monomers in a copolymer. We find, 



1 <«»(q)« ) (-q)>0- 1 '-'^ [ l J**M-W*n 



Jo \ s m ( 



V v ' a v ' p N 1-a Jo V sin(qlx/2) 

x cos(ql[(l - f)N + l]x/2) 



(1.2) 



T7 < 4>B (WB (-q) >0= —fj / dx . 2 - ; (L3) 

V Po^ Jo sin (qlx/2) 



y < A n) (q)^ } (-q) >o= ^ < ^ n) (q)C } (-q) >o= 
1 1 - a fN f 1 ( (l 3 2 \ sin(qr/(l - f)Nx/2) 



p N 1-a 



\ < 4"'(q)Qi?(-q) >o=^ < ^"'(qlQif'f-q) > = 



> N Jo 



(] ,l_S x 2\ sm 2 (ql(l-f)Nx/2) 



(1.4) 
(1.5) 



PoAJo ^4 4/ sin 2 (g/a;/2) 

i < Qif^Qifi-q) > = 1 < C } (q)C } (-q) >o= 

1 f\x(^(i-xr-l(i-^ + ^ sin2{ql{1 - f)Nxm 



(1.6) 



2p Ay V!6 2 V 7 2; sin 2 (g/x/2) 
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I<Q<f)(q)Q<f>(-q)>„= 



2,ooiV./o \!6 2 V 7 2; sin 2 (g/x/2) 

Finally in Eqs.(I.l-I.7) we take a limit of </ — > 0, N A ,N B — > 00 such that 
N A /(N A + N B ) =const and g/iV =const. 
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Figure Captions 

Fig.l The schematic picture of the unit cell of cubic bicontinuous double diamond 
phase formed in AB diblock copolymers. The junction of the A-B blocks are 
located at the surface shown in the figure. The surface defines the boundary be- 
tween A rich and B rich domains. Symmetry Pn3m. The experimental studies 
suggest that the surface describing the boundary of domains in diblock copoly- 
mers is of constant mean curvature and is determined by the area minimization 
subject to fixed volume fraction of A and B blocks. 

Fig. 2 The schematic picture of the unit cell of gyroid structure formed in AB diblock 
copolymer systems. The junction points of the A-B blocks are located at the 
surface shown in the figure. Symmetry Ia3d. 

Fig. 3 Phase diagram of the polystyrene-polyisoprene system 30 ' 31 ). The bcc phase is 
the body center cubic phase, where droplets of the minority component in the 
diblock copolymer form an ordered lattice. The hex phase is the hexagonal 
phase where the parallel cylinders of the minority component surrounded by 
the majority component form a hexagonal lattice. The obdd phase is the or- 
dered bicontinuous double diamond structure consisting of two channels each 
of diamond symmetry separated by the surface (Fig.l). ODT stands for order- 
disorder transition, x is the Flory interaction parameter, N is the polimeriza- 
tion index and / is the fraction of one component in the diblock copolymer 
system. The region of stability of the gyroid phase (Fig. 2) which has been 
recently observed in this system is not shown here. After Ref. 28. 

Fig. 4 Bridge and loop configurations in the lamellar phase of multiblock copolymer 
system. After Ref. 44. 

Fig. 5 (a) On the left is shown the permitted configuration of an unknotted ring. It 
cannot turn into the knotted configuration shown on the right (not vice versa). 
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(b) Two unknotted rings are shown in the allowed configuration. It cannot 
turn into the knotted configuration shown on the right. After Ref.46. 

Fig. 6 Microscopic model of an AB diblock copolymer molecule with different stiffness 
of the blocks. The size of the A segment, bA, is much smaller than the size of 
the B segment 65. After Ref.52. 

Fig. 7 The one- loop diagrams contributing to the two-body vertex function. The solid 
lines represent while the higher order vertex functions are represented 

by points in the diagrams. The n-body vertex function has n lines emanating 
from the point representing this function. The dashed lines are not parts of the 
diagrams but are included only for clarity (according to Ref.(7)). Apart from 
l/r^ , only and appear in the structure of the diagrams. 

Fig. 8 The one-loop diagrams contributing to the three-body vertex function. The 
legend as in Fig. 1. Apart from l/T^ only and appear in the 
structure of the diagrams. 

Fig. 9 The plot of the logarithm of the scattering intensity versus the logarithm of the 
temperature difference. After Ref.ll. Close to the critical point the intensity 
diverges according to the power law \T — T c |~ 7 . The experiment was performed 
here on the symmetric mixture of the PDMS and PEMS. 
Fig. 10 The inverse of the scattering intensity at zero q vector. The temperature T x 
marks the Ginzburg region where there crossover from the mean-field to the 
non-mean-field behavior. Tf, is the critical temperature on this plot. Please 
note that T Sjm / (mean-field temperature) is not correctly identified (see text 
and Ref.15,74). After Ref.10. The experimental system studied in this paper 
is PI-PEP mixture. 

Fig. 11 The scattering intensity 22 ) for the diblock copolymer in the random phase 
approximation for the system with composition / = 0.25 as a function of 
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q 2 R 2 , where R 2 = Nl 2 /6 is the radius of gyration; solid line (xN = 12.5), 
dashed line (xN = 16.0) and the dotted line (xN = 17.5). The peak in the 
scattering intensity is fixed and does not change with temperature, which is 
here proportional to ~ 1/x- 
Fig. 12 The phase diagram for the diblock copolymer system in the random phase ap- 
proximation (after Ref.22 and Ref.23). The diagram is symmetric with respect 
to f=0.5 (symmetric diblock copolymer). The transition to the body centered 
cubic (BCC) and hexagonal (HEX) structure is first order. The direct second 
order transition from the disordered phase (DIS) to the lamellar phase occurs 
only at / = 0.5. 

Fig. 13 The phase diagram for the diblock copolymer system in the self consistent 
one-loop approximation (Hartree) (after Ref.23) for N = 10 6 . Compare with 
Fig. 12; in the limit N — > oo the results of the Hartree approximation reduce 
to the random phase approximation. The dashed line marks the divergence of 
the scattering intensity in the random phase approximation (of course in the 
one-loop approximation for finite N the structure factor does not diverge at 
this line). 

Fig. 14 Peak position, q* , (in A -1 ) versus the polimerization index, N for the PEP- 
PEE diblock copolymer system at fixed temperature T=296 K. After Ref.20. 
The vertical solid line indicates the microphase separation obtained in the 
theory. The solid line of slope -1/2 indicates the result of the random phase 
approximation. Open squares are the predictions of the one-loop approximation 
theory and full triangles are the experimental results 78 ") . 

Fig. 15 The phase diagram obtained in the self consistent field theory. Here C is the 
cubic bcc phase, H is the hexagonal phase, L is the lamellar phase and G is the 
gyroid phase (Fig. 2). The diagram is symmetric with respect to reflection at 
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f=0.5. After Ref.24. The transition from the disordered to the lamellar phase 
is second order; exactly as in the random phase approximation. It is marked 
by a dot in the phase diagram. All the other transitions are first order. 

Fig. 16 The phase diagram of the diblock copolymer system with conformational asym- 
metry. The size of the A monomer, I a, is equal to y/TOls, where Ib is the size 
of the B monomer. The legend as in Fig. 15, where I a = Ib- Please note that 
conformational asymmetry does not change the topology of the phase diagram 
in this theory. After Ref[25]. 

Fig. 17 The scattering intensity, calculated in the random phase approximation, for 
the symmetric ring diblocks and linear diblock copolymer system at infinite 
temperature % = 0. The scattering from rings is less intense and the peak is 
at higher q vectors. After Ref.48. 

Fig. 18 The instability towards the ordered phase (divergence of the scattering inten- 
sity). The transition temperature (~ 1/x) is approximately 40% lower in the 
rings than in the linear diblocks. After Ref.48. 

Fig. 19 Comparison of the radii of gyration for linear and ring (cyclic) symmetric di- 
block copolymers of length N = 20. After Ref.47. Here T c is the temperature 
of the microphase separation (to the lamellar phase); e is the enegy interaction 
parameter between the monomers and a is the monomer size. The inverse of 
this temperature, in the random phase approximation, is shown (Fig. 18) for 
/ = 1/2- 

Fig. 20 The distance between the centers of masses of A and B blocks in the ring diblock 
copolymer. After Ref.47. It indicates that the rings strongly stretch as the 
temperature is lowered towards the temperature of the microphase separation, 
T 

Fig. 21 The radius of gyration for the single block in the ring copolymer a = A, B. The 
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blocks shrink as we lower the temperature , similarly as in the case of linear 
diblocks. T c is the temperature of the microphase separation. After Ref.47. 

Fig. 22 Phase diagram for the ABA triblock copolymer system calculated in the ran- 
dom phase approximation. Here BCC is the body center cubic phase, HCP is 
the hexagonal phase and LAM is the lamellar phase. After Ref.84. The two A 
blocks have the same number of A monomers. 

Fig. 23 Phase diagram of the tenary mixture at equal A and B homopolymer concen- 
trations. Here only instability lines are marked. The three phase coexistance 
region is marked only schematically. In the region named lamellar there should 
be more different ordered phases, but they have not been studied, here. Also 
note that the point where all lines meet is the Lifshitz tricritical point. After 
Ref.35. See also Ref.34 and 36. 

Fig. 24 The same phase diagram as shown in Fig. 23. Here the lines in the disordered 
phase shows the various stages of local ordering. LD denotes the disorder line, 
LF is the Lifshitz line for the correlation function and LE and EG are the 
Lifshitz and the equimaxima lines of the rjr] structure function. L is the Lifshitz 
tricritical point. The region between the LD and LF shows the region where 
the correlation function starts to oscillate showing a local ordering. After line 
LF the ordering becomes dominant and the peak in the structure functions 
appear. The region between the LG line and LE-EG is the region where the 
homopolymers order and where their structure functions develop a peak. At 
the EG line the homopolymer homopoymer structure function shows two peaks 
et equal height, one at zero q vector and one at non zero q vector. After Ref[35]. 

Fig. 25 The stability limits of the disordered isotropic phase (I) against lamellar (L) 
(solid line) or nematic (N) perturbations (dashed line) in the rigid-flexible di- 
block copolymer system for (a) vbb/x = 0, (b) vbb/x = 0-4, (c) vbb/x = 3-0. 



Dot denotes the trier itical point. After Ref.50. 

Fig. 26 The structure function G^q,(q z ) (case ((a),(c)) and its one dimensional Fourier 
transform (density-density correlation function) G^^(z) (case ((b),(d)). Here 
N is the number of monomers in a diblock copolymer, po is the monomer density 
and N B l = (l-f)Nl is the length of the rigid part. (a),(b)f=0.2, %iV =19.51, 
and vbb/x = (dashed line), and f=0.2, xN =19.51, vbb/x =0.4 (solid line). 
(c),(d) f=0.62, v B b/x = 0, %iV=11.75 (dashed line) and f=0.62, v BB / x =3.0, 
xN— 11.3b (solid line). The parameters for the solid line are close to those of 
the tricritical point shown in Fig. 25. After Ref.50. 

Fig. 27 The structure function Gqq((/ 2 ) versus q z and its one dimensional Fourier trans- 
form (nematic order parameter-nematic order parameter correlation function) 
Gqq(z) versus z for f=0.62. The parameters for the solid and dashed lines are 
as in Fig(26c,d). After Ref.50. 

Fig. 28 The structure function G^Q(q z ) versus q z and its one dimensional Fourier trans- 
form (density-nematic order parameter correlation function) G^q(z) versus z 
for the same two systems as in Fig. 26; (a),(b) f=0.2; (c),(d) f=0.62. After 
Ref.50. 
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